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Prefaţă 


Acest manual își propune să introducă pe elevi în metodele de calcul apro- 
ximativ, atît de utile astăzi. | 

Orice problemă practică, pentru a fi rezolvată pe calculator trebuie să fie 
adusă la un algoritm programabil și apoi algoritmul să lucreze. 

Ne-am străduit să punem în evidență pe tot parcursul acestui manual 
diferitele tipuri de erori ce apar în calcul: cu efectele lor asupra rezultatului. 

De fapt, unul din aspectele esenţiale ale acestui obiect este analiza erorilor. 

Acest lucru nu este deloc uşor ci dimpotrivă. 


În consecință $3 dir capitolul III va fi studiat facultativ întrucît are un 
grad de dificultate destul de mare. 


L-am inserat totuși în text pentru a sprijini punctul de vedere anterior 
enunțat. [3 : 

Peste tot am folosit interpretarea geometrică a faptelor în scopul înțelegerii 
prolunde (și mai ușor) a fenomenelor prezentate. 

Pentru ca munca celor ce folosesc acest manual să nu fie în zadar, va trebui 
să se efectueze multe exemple numerice. 

Aceste exerciții îi vor forma pe elevi în „spiritul“ „convergenţei numerice“. 
Vor fi în măsură să privească niște rezultate scrise de calculator și să-și dea 
seama dacă sînt bune, iar dacă nu, să descopere de ce nu sînt bune. 

În aplicațiile matematicii un anume „simț“ al numerelo:: este foarte util 
şi se poate cîştiga numai prin ducerea, la bun sfârșit a multor probleme de calcul. 

Orice sugestii din partea cititorilor vor fi primite cu interes. 


Dedicăm acest manual profesorilor noștri Paul Tăvăluc din Oradea și Roșu 
Dumitru din București. 
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Capitolul | 
ERORI. CONCEPTE DE BAZĂ ÎN CALCULUL NUMERIC 


$ 1. Introducere 


Calculul aproximativ joacă un rol central în matematica aplicată. 

În practică ne mulțumim foarte adesea cu o valoare aproximativă a can- 
tității pe „care o găsim“, alteori ignorăm unele cantități mici în comparaţie cu 
altele, iar nu de puţine ori „datele de intrare“, ale problemei puse, sînt inexacte, 

Introducem următoarele notații (nu sînt exacte din punct de vedere 
matematic) utile în practică: 

a & b (sau b > 4) citim: „a este mult mai mic decît 6“ (sau „b este mult 
mai mare decît 4). 

Evident, acest „mult mai mare“ sau „mult mai mic“ depinde de contextul 
problemei (uneori ab dacă a< - - 5, în timp ce în alte situații a&b 


dacă Pa IE i - 
zi 500 


az b; „a este aproximativ egal cu b“ și înseamnă că |ja—b|&c, 
unde c este ales de contextul problemei. 


$ 2. Surse de erori 


Orice rezultat al unui calcul numeric poate fi influenţat (uneori dezas- 
truos !) de diferite tipuri de erori. 

a) Erori în datele de intrare. De exemplu dacă coeficienţii unui sistem 
liniar sînt fracții ordinare, la transformarea lor în fracții zecimale se introduc 


erori [-: = 0,3333339), Alt exemplu, considerarea în calculul cu numere ira- 


ționale (m, 2) a unui număr finit de zecimale. 

Uneori datele de intrare sînt rezultatele unor măsurători carc sînt siste- 
matic afectate de erori. A 

În problemele tehnico-științifice (rezolvarea cu calculatorul a ecuaţiilor 
diferenţiale, de exemplu) apar foarte des sisteme liniare algebrice în care 


! 5 


|! 


coeficienţii sînt! integrale ce trebuie calculate aproximativ şi deci coeficienţii 
sistemului liniar sînt afectaţi de erori. 


b) Erori de rotunjire în timpul calculului. Dacă facem un calcul utilizînd 


numai cifre în toate numerele cu care lucrăm, atunci produsul exact a două - 


numere este format cu 2/ sau 2/ — 1 cifre. Dar noi nu putem utiliza la pasul 
următor de calcul decît rotunjirea la 7 cifre. Astfel se produc erori care pe par- 
cursul unui lung proces de calcul pot afecta serios rezultatele. 

c) Erori de trunchieve. Acestea sînt erori ce provin din limitarea la o „parte“ 
din termenii unui proces de calcul înainte de a. ajunge la limită. De exemplu, 


: a PS ara l | | 

din sumă progresie! geometrice infinite 1 e pe e SRR 
Di 3 22 198 

numai primii patru termeni, diferenţa dintre suma progresici și acest rezultat 

se numeşte eroare de trunchiere. 


reținînd 


Tot o eroare de trunchiere facem şi atunci cînd aproximăm 'derivata unei 
funcţii f(%) într-un punct x prin raportul: 


unde /1 este „mic“. 


De asemenea același tip de eroare se face cînd pe o porțiune mică în jurul 
unui punct aproximăm o funcție cu tangenta la graficul său în punctul res- 
. 
pectiv. 
d) Erori de simplificare a modelului şi erori comise de om şi mașină. 


Uneori, în practică, în unele modele, idealizăm anumiți parametri pentru 
a putea prelucra modelul. 


De asemenea se poate greşi în programul pentru calculator sau în perfo- 
rarea cartelelor, erori de operare etc, 


$ 3. Erori relative şi erori absolute 


Fie ă o valoare aproximativă pentru cantitatea a (a este valoarea exactă) 
avem: 


Eroarea absolută în ă este 3 —a. 
Eroarea relativă în ă este (ă — a)/a daci a 20. 


Foarte des eroarea relativă se dă în procente. De exemplu, 4%, înseamnă 
o eroare relativă de 0,04. 


Acum „a — ă“ este corecţia ce trebuie adăugată lui ă pentru a obțin 
valoarea exactă, : 

Corecţia și eroarea au aceeași valoare absolută dar semne contrarii. 

Eroarea poate fi pozitivă sau negativă dar noi putem introduce o mar- 
gine pentru eroare care este întotdeauna un număr pozitiv. 

Dacă a = ă 4 e înseamnă că |ă —a| se 

De exemplu a == 0,462 E 0,011 înseamnă: 


0,451 s a s< 0,4% 


$ 4. Sisteme de numerație. Reprezentarea numerelor 
în virgulă mobilă 


p zi SA b ae A 
În general sîntem obișnuiți să gîndim și să calculăm cu numere în baza 10. 
Astfel numărul 257,424 reprezintă: 


2: 1024 5 10147 10044. 101 | 2 + 102 4-4 - 103 


Orice număr real în baza 10 are o reprezentare unică, în forma de mai 
sus. Se exceptează cazul în care, în scrierea pozițională a unui număr zecimal 
(în baza 10) cifra 9 apare de o infinitate de ori. 

Exemple: 

1) numărul 99,9999 ... peprezintă același număr ca 100 

2) numărul 2,219999 ... reprezintă același număr ca 2,22, 

În general pentru reprezentarea unui număr real putem utiliza o bază 
oarecare B (număr natural) cu B>2 și atunci se poate arăta (cu suecpiii ca 
mai sus) că orice ntimăr real pozitiv admite o reprezentare unică de forma: 


=] =) 
Ap Br apa BO a: Ba hau: Birap- B+... 


unde coeficienții „a“ se numesc cifre în sistemul cu baza B şi sînt numere 
întregi pozitive cu proprietatea că: ' 


„ 0Osa,sB-—l. 


Avantajul mare al scrierii unui număr într-o bază care conduce la o scriere 
pozițională de forma: 


Adlmă «ee didos dq goe 


unde 4,4 1 +:- do reprezintă partea întreagă iar O, a ap... partea frachionară 
a numărului real, este că putem da reguli simple pentru operațiile aritmetice. 
Cu cît baza este un număr natural mai mic, cu atît aceste reguli devin 
mai simple. Acest argument stă la baza faptului că cele mai multe calcula- 
toare operează în baza 2 (în sistem binar cu două cifre 0 și 1). 
Amintim regulile de adunare»și înmulțire în sistemul binar: 
! 


0+0—0 0-0—0 
(0 il 03 0, 
1+0—1 1-0=—0 
jei2l 210 fala li 


i i ilizate sînt cele zu cifre de la 0 la 7 

Alice sisteme de numerație utilizate sint cele octale cu cilre de | 

și hexadecimale cu cifre de la 0 la 9 și literele A, B, C, D, E, [! de la zece 
la cincisprezece. 


De exemplu: . 


(13,25) = (101,012 = (15,2)e = (P,4)as 


A 
unde prin (no Oa (a (he am notat baza de lucru. 


Virgulele utilizate aici separă partea întreagă de partea fracționară ȘI 
se numesc respectiv: virgulă zecimală, virgulă binară, virgulă octală, virgulă 
hexadecimală. 

Un calculator este echipat de regulă cu două tipuri de operaţii aritme- 
tice, anume calculul cu virgulă fixă și calculul cu virgulă mobilă.. Virgula 
este zecimală dacă baza este 70, binară dacă baza este 2 etc. 

În cele ce urmează ne vom ocupa în special de baza zecimală. 

Calculul cu virgula mobilă zecimală se efectuează cu un număr constant 
de cifre după cum se va observa în cele ce urmează, 


Eat Pl a ut 
Ne oprim acum mai îndeaproape la” calculul cu numere reprezentate în 
virgulă, mobilă. 


Numim reprezentare zecimală normalizată cu virgulă mobilă a numărului 
A o reprezentare de forma: 


a = m 10, unde 0,1 s |mj<l 


£. 
m este număr real pozitiv sau negativ, 
9 este număr întreg pozitiv sau negativ. 


O astfel de reprezentare este totdeauna posibilă, pentru 440. 

Variabila m poartă numele de mantisă şi este o parte fracționară 
se numește exponent și este întreg. 

De exemplu numărul 542,36 va fi reprezentat ca: 

0,54236 - 105, deci 0.54326 este mantisa iar 3 este exponentul. 

Alte exemple: 

1) 4527 = 0,4527 - 104 mantisa este 0,4527 și exponentul 4 

2) 0,0002617 — 0,2617 . 10-2 mantisa este 0,2617 și exponentul —3 

3) — 42,15 = — 0,4215 - 102 mantisa este — 0,4215 şi exponentul 2 


„iar q 


4) — 0,012 = — 0,12 - 10”! mantisa este — 0,12 și exponentul —1 


Un număr zecimal în virgulă mobilă este memorat, în multe calculatoare, * 
cu mantisa și exponentul în aceeași celulă de memorie. 

Celula de memorie însă nu are prevăzută decît o singură poziție binară 
pentru semn și astfel se ridică problema memorării cînd și mantisa și expo- 
nentul sînt negativi. Soluția este în a aduna exponentului o constantă pozi- 
tivă, de care se ține seamă, în calculele aritmetice, 

De exemplu, într-o mașină zecimală cu două cafre  Bentru exponent vom 
aduna 50 la toți exponenții (ceea ce înseamnă căi exponenții din exemplele de ma: 
Sus vor fi reprezentaţi în interiorul calculatorului ca 54, 47, 52, 49), 

pe două poziții zecimale rezervate exponentului putem memora 100 numere 
diferite (numerele de la 0 1a 99). Cum exponentul poate lua și valori negative 
vom considera un interval de valori de lungime 100 simetric în jurul orig 


50, +49] 


inii: 


Deci exponentul real va lua valori în intervalul [— 50, +49] însă 
calculator va fi memorat (pe cele două poziţii rezervate lui) un număr cu 
mai mare. 

Într-un calculator numărul de cifre pentru g și m este limitat. 


Ca > 
o 


Să studiem acum operațiile în virgulă mobilă. 

|) Să presupunem că dorim să adunăm numerele: 

34 = 112,5 în virgulă mobilă memorat ca 0,1125 - 109 

Ya = 15 în virgulă mobilă memorat ca 0,15 - 102 

Așa cum sînt reprezentate numerele nu putem aduna mantisele deoarece 
virgula numerelor nu este aliniată. 

Diferenţa dintre cei doi exponenţi este |, vom deplasa atunci mantisa 
celui de al doilea număr (numărul cu exponentul mai mic) spre dreapta cu 
un loc, adică scriem: 

o = 0,015 - 108 


și acum putem aduna: 
3 da = (0,1125-4- 0,015) : 10% = 0,1275 - 105 


Observăm că am lucrat aici cu o mantisă formată cu 4 cifre. 
2) Să considerăm un alt exemplu: 


d -F Xa — 0,1213 - 1014 0,4691. 10-1, 


Deplasăm din nou mantisa numărului cu exponent mai mic spre dreapta 
cu un număr de poziții egal cu diferența exponenților, deci în acest caz cu 2 
ȘI obținem: 
Xu-F Xa = 0,1213 - 1014+0,004691 - 101 = 0,1259 . 101. 


Se observă că am pierdut cifrele 9 și 1 ale mantisei deplasate în afara 
capacității de calcul (deoarece am considerat mantisa ca avînd patru cifre). 
3) Presupunem acum că avem de efectuat următoarea adunare: 


0,9917 - 109+ 0,1230 . 105. 


Numerele de adunat au același ordin de mărime (același exponent) deci 
singura operaţie ce urmează a se efectua este adunarea mantiselor: 


0,9917 320, 1230) 251, 1475 


Rezultatul adunării mantiselor este un număr mai mare ca 1. 
Noi am definit însă mantisa ca o fracție zecimală subunitară, 


0, < |pl|<l. 


Operația care urmează a se efectua este deplasarea mantisei spre, dreapta 
cu o poziție zecimală și adunarea la exponent a unei unități. Această operație 
în . 
se numește renormalizare. 


Deci: Ș 
0,9917 - 105--0,1230 : 10% — 0,1114 - 107 


(din nou prin deplasare spre dreapta am pierdut cifra 7 deoarece calculul 
se efectuează cu 4 zecimale). 
4) Fie acum de ejectual produsul a două vumere Xa ȘI Xa în dDirgulă mobilă: 


0,1763 - 102 x 0,2036 . 105. 
Normal, se înmulțesc mantisele, 


Ma * Ma — 0,1763 x 0,2036 = 0,0358, 


iar exponenţii se adună qi + ga = 3. 
Deci, 


A+ Xa = 0,0358 - 108, 


Nici de data aceasta mantisa (în modul) nu este cuprinsă în intervalul 
[0,1; 1). 


Operația de renormalizare constă în acest caz în deplasarea spre stînga 
a mantisei 0,0358 cu o poziţie zecimală și scăderea unei unități din exponent 
(adică aş * 42 — 0,3580 - 102), 


Observăm deci că. operaţia de renormalizare, în urma efectuării unci ope- 
rații aritmetice constă în „aducerea“ mantisei în intervalul [0,1; 1) cu „ajus- 
tarea“ corespunzătoare a exponentului. a 


Observatie. În majoritatea calculatoarelor înmulțirea se efee- 
tuează pe registru de memorie dublu (adică o celulă dublă de 
memorie). Dacă se înmulțesc două numere de lungime „!*, rezul- 
tatul de „24* sau „24 — 1“ cifre este obținut în acest registru de 
lungime dublă. Dacă după efectuarea înmulțirii este necesară 
venormalizarea. (numai cu deplasarea spre stînga în acest caz) 
cifra de pe ultima poziţie (poziția „/*) a numărului va proveni din 
a doua parte a registrului dublu de memorie şi va fi deci corectă. 
În: exemplul de măi sus 44 + Aa = 0,03588468 - 102 în registrul - 
de lungime dublă iar rezultatul renormalizat pe patru zecimale 
va fi 0,3588 . 102, 

Dăm în figura 1.1 organigrama „ădunării a două numere în virgulă mo- 
bilă, : 

Amintim că, în organigrama din figura 1, un dreptunghi reprezintă un bloc 
de calcul și de atribuire (se efectuează operaţiile din membrul drept al sem- 
nului egal și valoarea obţinută se atribuie parametrului din partea stingă 
a semnului egal) iar un romb reprezintă un bloc de comparăție și de decizie. 


Presupunem acum că avem un calculator în care sînt prevăzute 5 poziţii 
zecimale pentru mantisă și două poziţii zecimale pentru exponent, 


Cel mai mare număr normalizat ce poate fi memorat cu aceste facilități 
va fi: ; 


9,99999 - 10% în realitate 0,99999 . 104 deoarece am convenit să adunăm 
exponentului real numărul 50. 


Iar cel mai mic număr normalizat (în modul) va fi: 


0,10000 - 10? în realitate 0,10000 - 10-50 


Normalizează, m3 
m =m3 1043 


Fig. 1. 1. Orgauigrama adunării a două numere în virgulă mobilă. 


(3 = Ma BO, Xa — Ma Be, x Xa — ma Br), 


Cînd printr-o operație aritmetică între doi operanzi am depășit una sau 
alta din limitele arătate mai sus (un număr mai mare de 1,0: 10% — (depă- 
şire ilotantă superioară) sau un număr mai mie de 1,0 - 10-% — (depășire 
flotantă interioară)) calculul se oprește. 


Exemple: 
1) 0,12641 - 10-45 x 0,98763 + 10-24 
'2) 0,45631 - 10%2 x 0,18764 . 1015 


vor cauza oprirea procesului de calcul. 


$ 5. Rotunjirea şi roturijirea prin tăiere (trunchiere) 


Există două căi de rotunjire a numerelor la un număr de zecimale („£*) 
Hxat anume, rotunjirea prin eliminarea tuturor zecimalelor care se găsesc 
la dreapta poziţiei £ (rotunjire prin tăiere) și a doua cale este rotunjirea obiş- 
suită adică printre numerele ce se pot forma cu „t“ zecimale alegem numărul 
cel mai apropiat de numărul dat. 


Exemple 


Rotunjirea şi votunjirea prin tăiere (trunchiere ) la trei zecimale în numerele 
de mai jos se face după cum urmează: 


1) 0,1597 se rotunjește la 0,160; 
2) 0,2456 se rotunjește la 0,246; 
3) 0,1574 se votunjeşte la 0,157; 
4) 0,23652 se rotumjește la 0,237; 


se rotunjeşte prin tăiere la 0,159 
se votunjește prin tăiere la 0,245 
se rotunjeşte prin tăiere la 0,157 

se votunjește prin tăiere la 0,236, 

Un mare număr de calculatoare utilizează rotunjirea prin tăiere în rezul- 
tatele obținute după fiecare operaţie aritmetică, însă dat fiind faptul că numă- 
rul de cifre utilizate în calcul este mult mai mare decît numărul de cifre în da- 
tele problemei nu se ivesc neplăceri. 

Aflarea marginii erorii relative pentru un număr aproximat prin tăiere. 
Să presupunem că avem la dispoziţie un calculator cu 4 cifre disponibile pen= 
tru mantisă. 

Să adunăm Xa — 0,1142 - 101 cu a =— 0,2190: 101, 

Obținem: 

13 — 31 -F dtp — 0,1142 - 1024 0,002190 - 101 == 0,11639 - 101, 

Acest număr se poate exprima sub forma: 

aa =—0, 1163-10 70,9. 108. 

Exponentul celui de al doilea termen este mai mic cu 4 decît exponen- 
tul primului termen. Întotdeauna cînd lucrăm cu o mașină cu patru cifre în 
mantisă, exponentul celui de al doilea termen trebuie să fie mai mic cu patru 
decît exponentul primului termen. i 

Deci dacă avem un calculator cu „“ cifre în mantisa unui număr repre- 
zentat în virgulă mobilă atunci orice rezultat sal unui calcul (obținut prin 
aplicarea oricărei din cele patru operaţii aritmetice) se poate reprezenta (evi- 
dent înainte de rotunjire) ca: 

b = my: NOT my, » 10%% 


unde m, este un număr format cu „i“ cifre. 
Valorile posibile pe care le poate lua 7, sînt: 


(al ese jegos 
în timp ce: 
(0) ese) e île 
Spunem că facem o aproximaţie prin tăiere (trunchiere) dacă, pur și 
simplu nu-l luăm în considerare pe m, și acesta este cazul practic cel mai des 
întîlnit la calculatoare. 
Eroavea velativă maximă se obiine cînd my este cel mai mare şi my cel mat 
mic. 


Valoarea maximă posibilă a lui my este mai mică decît | sar valoarea mi 
nimă a lui m, este 0,1. 
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Deci limita maximă a erorii relative va fi: 


1710 er 
0,1 - 10% 
La rotunjire; prin acelaşi procedeu, obținem o eroare relativă dată de 
5. ge E 
0,5100 = 0,5: 10-t. 
0, - 10% N 


" În general pentru o bază arbitrară B avem adevărat următorul criteriu: 
Presupunem că lucrăm în baza B cu „t“ cifre în mantisă (cifra binară 
care ne dă semnul numărului nu este socotită), Atunci orice mumăr real în 
vivgulă mobilă poate fi veprezentat cu o eroare velativă care nu depășește unita- 


tea mașinii definită de: 


Încă a =. 
: le, BI dacă utilizăm rotunjirea 
14. =— 


lu dacă utilizăm tăierea. 


În cele mai multe calculatoare „4“ se află între 10-€ și 10-15, 

Să presupunem acum că operanzii cu care lucrăm ar fi exact reprezen- 
taţi în mașină. Nu este sigur că rezultatul va [i exact reprezentat (să ne gîn- 
dim doar la faptul că rezultatul înmulţirii a două numere de lungime Mare 
„24 cifre). N 

Vom nota cu fi(ă +-y), f(x —y), Jilx-y), fllxly) rezultatele opera- 
țiilor de adunare, scădere, înmulțire și împărțire în virgulă mobilă pe care 
maşina le memorează în urma efectuării acestor operații și aproximării prin 
tăiere sau rotunjire. y 

Conform criteriului precedent asupra erorii relative avem 


a ON= (00% 
(309) 
unde (2) este una din operațiile +, —,.,/ 
Să răspundem acum la întrebarea: 4) Care este wmilatea maşinii pentru 
ui calcul cu virgulă mobilă cu 5 cifre zecimale pentru mantisă ? 


e 0,5 - 104 dacă se utilizează rotunjirea 
2 1 - 10-4 dacă se utilizează tăierea. | 


b) Fie acum x, = 0,12347 - 10% un număr rezultat din calcul. 

Case este marginea penlru eroarea absolută în acest calcul dacă s-a utilizat 
rotunjirea obișnuită? 

Eroarea absolută e, este mărginită de e, = 1 : 10%:0,5: 10%=—0,5.- 10523 

c) Pentru mumerele rezultate din calcul Xa = 0,12347 - 101% și axe=— 
— 020006 : 1010 erorile absolute vor fi mărginite de: ex — 0,5 - 101 şi ex = 
= 0,5 - 10, 


“3 
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d) Care este marginea pentru eroarea absolută în calculul celor trei mumere 
. - î: . .. . . >] 
de mai sus dacă s-a folosit rotunjirea prin tăiere : 


e =0,L-101 

ep = 0,1 + 1014 € 

eg = Oul - 10% 

Putem cere variabilei x, o precizie absolută de 0,1 : 10-52 Dar o precizie 
de, 0,1 - 10-40? Evident că nu, deoarece valorile de 0.1 : 10-6 și O, + 10-40 


sînt mult prea mici în comparație cu eroarea cu care calculatorul ne „furni- 
x x 2 
zează“ acest număr (0,5 » 10%. 


Observaţie. La efectuarea unui calcul cu virgulă mobilă trebuie 
să avem permanent în vedere că ultima cifră zecimală a mantisei 
este afectată de erori de rotunjire sau trunchiere. Deci dacă am 
obținut dintr-un calcul cu virgula mobilă, cu 5 cifre zecimale 
pentru mantisă, numărul 4 = 0,23483 : 107 marginea pentru eroa- 
rea absolută este de ordinul 0,000005 - 107 — 0,5: 102 în cazul 
cînd se aplică rotunjirea sau 0,00001 : 107 = 0,1: 10% cînd se 
aplică: truchierea și deci nu putem cere o precizie mai bună decît 
0,5 x 102 respectiv 0,1: 105. 


Acum vom arăta că asociativitatea în adunarea în virgulă mobilă nu se 
păstrează. 


Presupunem că lucrăm cu o maşină cu şapte cifre în mantisă. Fie: 
a = 0,1234567 109 i 
= 0,4711325 + 104 


= 


c=—b 
avem: 
JI + c) = 0, Ja + /1(6+ c)) = /l(a) = 0,1234567 - 109 
pe de o parte, iar pe de alta: 
Jl(a + 5) = 0,00001234567 - 1044 0,4711325,- 101 = 0,4711448 : 101 
fi(fi(a-+ b) + c) = (0,4711448 — 0,471 1325))10% = 0,0000123 - 10% = 
= 0,1230000 - 100, 


$ 6. Propagarea erorilor 


Dacă 4, = 2,32 + 0,03 ȘI Xa = 1,24 4+ 0,02 atunci care este marginea 
erorii pentru x — 9? 

Cea mai mare valoare posibilă a lui x, este 2,35 iar cea mai mică valoare 
posibilă a lui +2 este 1,22, deci cea mai mare valoare posibilă a lui 44 — Xa 
este 2,35 —'1,22 = 1,13. 

Similar cea mai mică valoare a diferenţei A — Xp este 1,03 
Deci: 
1,03 < 4, — 2 s 1,13 
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sau 


În cazul general dacă x, = 
Ă 
avem: 


pai ză [2 Wiraa (Xa Eu £) S XX Ss Sa Es (Ga == c2), 


Ș — Xa — (e. -F £2) S< Ma — Wa < 1 — Xa (21 + £2), 
( 


Prin inducție ajungem la următoarea: - NI 

Regulă. La adunare și scădere, marginile pentru eroarea absolută sînt 
date de suma marginilor erorilor absolute în operanzi. 

Conform definiţiei erorii relative (vezi $ 3) avem: 


a Pra Pi A, 
E Ea. ” (unde v este eroarea relativă) 


3 « 
lee 35 = AURA) 
Dacă î şi Xa au erorile relative 7, și 72 avem: 
a = Xa (L-A ra) all 2) = raa( LA ra) (Lara). 


Urmează că eroarea relativă în produsul 34, este dată, de relația: 


Xa — Xa 


Xa 


= AIE 0) (ET) A Aaa 


Dacă |ri| <£1 şi [ra] < | atunci. putem neglija termenul rara și obținem 
eroarea relativă în Xa ca fiind aproximativ 71 + 7a- 
De exemplu dacă 7 = 0,03 şi ra = — 0,02, eroarea relativă în produs este: 


PF pt ia = 0,03 — 0,02 — 0,0006 — 0,0094 x 0,01: 


iti 2 a e Spa a 
Să procedăm la fel pentru a găsi eroarea relativă la împărțire. 
"Avem: 


deci: 


de unde: 


Baa ați le te za 
Za Ya L-a 
iar 
Ei F 
pa a _l+r LL VA) 
Star A 1 7a 
Se 13 4 
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. 
Dacă |n] & | și |ra] < 1 obţinem: 
N — Va 
E OI ALA A A 
Lp 
Dacă marginile erorilor relative în 4 Xa Sînt pu și ea atunci put ge este 
cea mai bună margine superioară pentru |r+ ral ca şi pentru |n — ra] 
În concluzie avem următoarea regulă: 
La înmulțire și la împărțire marginile pentru erorile relative în operanzi 
se adună. 
O situaţie des întilnită în practică, care conduce la rezultate eronate este 
următoarea: 
Să Presupunem că avem, de făcul o scădere între doi operanzi, în care dife- 
vența este mult mai macă deci fiecare dintre ei. 
Fie y = i — Xa, Să notăm erorile în Xy ȘI Xe cu Axa Și Axe atunci: 


Ay] < |Az1+ Asa] 


deci 


Asi î IAzal + Idel 
y la — al 
| 


Aplicaţie: Dacă a = 0,6724 + — 104 şi 
2 


Na = 0,6723 2 — "104 
2 


atunci: 
44 — Xa = 0,0001 4 0,0001. 
Această anulare a rezultatului poate fi ocolită printr-o scriere convena- 
bilă a formulelor de calcul sau aplicarea unui alt algoritm. 


Să considerăm următorul exemplu: 
ps Lisi 
ecuația x — 56x+ | = 0 are rădăcinile: 


x = 28 — 7832 28 — 27,982 — 0,018 4 LL 10-3, 
2 


Se observă că deși rădăcina pătrată s-a calculat cu cinci cifre semnifica- 
A a : ete ab - ma E ce PA aa 
tive, totuși în prima rădăcină (4) primim numai două cifre semnificative iar 
în 4 eroarea relativă este mai mică decît 10-5. 

Cum produsul rădăcinilor este vaz = 1, putem utiliza următoarea cale 
pentru calculul lui a: 


= 0,01786288 


X= = - 
ap. 55,982 


cu o eroare relativă mai mică decît 105. 
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Așadar x = 0,0178629 + 0,0000002 și astfel cu aceeaşi estimare a rădă- 
cinii pătrate obținem pentru xy'cinci cifre semnificative (în loc de două) uti- 
lizînd informația dată de produsul rădăcinilor unei ecuaţii de gradul al doilea. 

În general dacă |3| < x atunci vom scrie: 
vă — 8 


Va 5 Va i: Ver 5+ya d 


De asemenea în multe situații este preferabil să utilizăm pentru calculul 
diferenței: 


cos (41-38) — cos x 


produsul, 
— 2 sin îi sin (+ + =) 
A, Z 


Dăm acum cîteva reguli simple în vederea executării unui calcul practic 
mai precis. 

I) Cînd se adună sau se scad numere (un şir de numere) să se înceapă 
cu cele mai Mici. “ 

II) Să se evite pe cît posibil cu putință scăderea a două numere pozitive: 
aproximativ egale. Deseori, după cum am văzut, putem rescrie formula de calcul 
pentru a evita acest lucru. 

III) O expresie de forma a(b — c) poate fi scrisă ab — ac iar (4—b)/c 
ca a/c — bjc. Dacă numerele sînt aproximativ egale să se facă scăderea înainte: 
de înmulţire evitînd astfel erori suplimentare de rotunjire. 

1V) Pe cît posibil să se minimizeze numărul operaţiilor aritmetice. 


E xercihia, 


7. Scrieţi în baza 10 numerele (1000), (100000). (464)s. 


2, Fie numărul 1/3 de memorat într-o celulă de memorie cu 4 poziţii 
zecimale. Cum apare acest număr în calculator? Cu ce eroare? Care este 
corecția? Același lucru pentru numerele 7 și i 25 

3. În, cazul reprezentării în calculator a numărului 1/3 (indiferent pe 
cîte poziții zecimale (în număr finit)), coincide eroarea de rotunjire cu eroarea, 
de trunchiere ? 


4. Pentru reprezentarea (cu rotunjire prin tăiere) a unui număr irațional! 


în calculator, coincide eroarea de rotunjire cu cea de trunchiere? 

5, Pentru un calculator ce lucrează cu 7 poziţii zecimale, care este 
unitatea ? 

6. Un calculator are rezervate 2 poziții zecimale pentru exponent. De 
ce se adaugă, în reprezentarea în calculator, numărul 50 exponentului real? 
Cum vor fi reprezentaţi pe aceste două poziţii exponenții —7, 6, —22, 11,20? 
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7. Avem un calculator cu 8 poziţii binare pentru exponent. Care este 
cel mai mare număr întreg ce „încape“ în cele 8 poziţii binare? Care este 
jumătatea acestui număr? Care este numărul ce se va adăuga exponentului 
real în reprezentarea în virgulă mobilă în calculator ? În ce interval se situează 
exponentul real? 


8. Fie de efectuat următoarele operații: 

a) 2,1330 - 7; b) 0,2567 -e;c) m - e; d) 0,2135—0,212; e) 1,2137 +82,2137 
pe un calculator cu două. poziții zecimale pentru exponent și 4 poziţii zeci- 
male pentru mantisă. : 

— Scrieți numerele în reprezentarea virgulă mobilă (normalizat) cu 
rotunjire. 

— Efectuaţi operaţiile pe celula simplă de memorie. 

— Renormalizați. 

— Efectuaţi calculele pe celulă dublă de memorie. 

— Renormalizaţi. ţ A 


9. Determinaţi eroarea relativă maximă pentru x, unde: 
e 2), 3220) 2 
Xp = 22,3 + 0,2, 
sii 203,3: 042: 


Care din aceste trei numere are precizia relativă cea mai bună? 
10. Calculaţi eroarea relativă maximă pentru: 
BE aa 
3 
unde xy, Xa, 43 sînt cei de la punctul precedent. 
1]. Să se efectueze, cu precizia cea mai bună, adunarea: 
Sat ah Xa 
în virgulă mobilă cu. $ zecimale (exponent 2 poziţii zecimale) pentru 
se AA des s 
Xa — 0,000043888, 
sa —.0,0082798. 


12. Calculaţi: 
1 ah 
? =j 1 iata 
- A Late) 
pentru n = 0, 1, 2, .... 8, utilizînd formula de recurenţă: 


1 
ala Î9a-a il 
VĂ 
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a) Cum rezultă formula de recurenţă ? 


b) Utilizaţi trei zecimale în calcul și rotunjiți. 
c) Arătaţi că ya > va şi ya <0 ceea ce este evident absurd, Explicaţi 


de ce? 


73. Utilizaţi în exemplul precedent formula: 


Dea Sa e 


Considerînd că vo & yp arătaţi că ys 2 0,019 și calculaţi vw, 63 sees o 


| 


Capitolul Il 
REZOLVAREA SISTEMELOR DE ECUAȚII ALGEBRICE LINIARE 


$ 1. Punerea problemei. Condiţionarea sistemelor 


Să considerăm sistemul de două ecuaţii liniare algebrice cu două necu- 
noscute: 
At Aa39 = ba (0) 
da + aa = ba. 


Coeficienţii necunoscutelor și, termenii liberi. îi considerăm numere reale. 
Vom nota cu A matricea, coeficienţilor necunoscutelor, deci: 


A Au da 


doi aa 


și să presupunem că determinantul matricei A (îl notăm cu det A) este diferit 
de zero, deci sistemul (1) admite soluţie unică. 

În această situație matricea A admite o inversă pe care o vom nota 
cu Al. 


Avem: 
aa — au] 
det 4  detA]| 
A = 
— Mai Au 
det A det 
ară 
Map — ial] 
a det A det A 
4-4 = âa||.| = 
da asa] | 2 aa 
|det A det 4 
| udaa —— dada 0 


i] 


det A | | | 


— Ay Fool 
det A 
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Deci A-A-1=I unde / este matricea unitate (la fel se verifică 
AF e Ale) 

Definiţie. — Se numește lungimea euclidiană a matricei A (o vom nota 
cu ||Allz), numărul pozitiv 


JA 


Calculînd ||A-1|y obţinem: 


ap, 


p = (aul! | aa2|2-+ |â2.]%A | 22|?) 


At 2 (zau PE Lâna? aa Pf aa PE 
| Ci |det A | 
Deci: 
în da [24 | 32 [24 | azi 24 laz2]* 
HA « at, = Iau E Lea E lea (2) 


ldetA| 


Definiţie. — Numărul pozitiv dat de relația (2) se numește numărul de 
condiționare euclidiană a sistemului (1) (îl vom nota cond, (4)= 
=> A le Ale). 


Observaţie. Există şi alte numere de condiționare pentru sis- 
teme liniare. Despre unele vom vorbi în acest context iar despre 
altele vor afla cei ce vor intra în studiul aprofundat al tehnicii 
de calcul. 


Din punct de vedere teoretic sistemul (1) admite soluție unică dacă. 
det A 4 0. Din punct de vedere al calculului practic lucrurile nu stau 
deloc astfel. 

Să considerăm următorul sistem de ecuații, pe care ne propunem să-l vezol- 
văm considevînd coeficienții săi dați cu două cifre deci, 


34 9y= 17 (3) 
js IEI = albe) 


Cu regula lui Cramer rezultă imediat că unica soluție a sistemului este: 
lips 

Din punct de vedere geometric, a rezolva sistemul (3) revine la a intersec- 
ta două drepte (fig. I[.1). 

Considerînd x» = 2,121, y = 0, în (3) obținem: 


[esa 9y = 16,968 


() 
9x-+-10y = 19,089. 

Prin rotunjire, la două cifre semnificative, membrul drept al sistemului 
(4) corespunde cu membrul drept al sistemului (3). În concluzie, soluția 
x — 2,121, y — 0 trebuie acceptată ca fiind la fel de bună ca cea unică 
întrucît coeficieriții și termenii liberi ai sistemului inițial au fost dați cu două 
cifre. 

Din figura dată se vede că dreptele sînt aproape confundate, iar punc- 


tul x — 2,121, y = 0 nu se găseste pe nici una din cele două drepte, dar el se 


găseşte foarte aproape de amîndouă. 
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x(2,121:0) i 


Fig. II. 1 


Calculind numărul de condiționare al matricei asocihte sistemului (3) 
obținem: 


condy(4) = dart Dn AU 326. 


A După cum observărn acest număr este destul de „mare“ în sensul următor: 
dacă efectuăm calculul în rezolvarea sistemului (4) cu trei zecimale exacte 


atunci obținem soluția + = 2,121, y = 0, pe cînd rotunjind în (4)'la două 
cifre avem soluția: 3 = |, y=1. 
Așadar numărul de condiţionare ne dă o oarecare imagine despre sensi- 


bilitatea soluției unui sistem liniar algebric la perturbații mici în membrul 
drept sau eventual în coeficienţi. 


În fapt, problema rezolvării sistemului. (4) poate fi privită și în telul 
„următor: la o perturbare foarte mică a membrului drept, vectorul 
| 17) | | 
| 19,098 
în soluție avem o perturbare destul de mare căci vectorul 
TI 121|) 
li od 


Să introducem acum pentru sistemul (1) și un alt număr de condiționare, 
Anume asociem matricei A, numărul pozitiv: 


devenind 


| 
a devenit | E 


IA] = max (| aaa | + laz2]; azil --|azal) 


(cu alte cuvinte valoarea maximă dintre numerele pozitive obținute prin 
adunarea pe linii a modulelor elementelor matricei A) | 


. 


Observaţie. — Se remarcă faptul că și acest număr pozitiv 
detinește o lungime pentru matricea A. 


; 

Ă De asemenea asociind și matricei inverse, după aceeași regulă, numărul 
pozitiv ||A4-1]|„, numărul de condiționare a matricei A (pe care-l vom nota 

: cu cond(4)) este prin definiție produsul numerelor pozitive ||4!|, și ||A-1l 
deci: Ă 


| Ala = ||AT||i = conda (4): 
De exemplu pentru sistemul (3) avem: 
Al = max (174 19)) = 19. 


0 
9 Sa 


Cum A-L= | sal avem! 


10|4+19]; )91+1—8]=19. 


[43]. = max (|: 
Rezultă că: 
erojara (04) = Ie) zl) sall 


Observăm că și acest număr de condiționare este de același ordin de mă- 
rime ca și cond, (4). 

Metodele de rezolvare a sistemelor liniare algebrice pe care le dezvoltăm 
| aici se justifică prin faptul că metoda lui Cramer de exemplu este foarte 
neeconomicoasă (calculul unui determinant necesită foarte multe operaţii) 
și practic utilizabilă numai pentru sisteme de două sau trei ecuaţii cu tot atitea 
necunoscute. 
| Problemele practice dau naștere la sisteme de ecuaţii de ordinul sutelor 

sau chiar miilor de necunoscute și atunci au trebuit căutate metode eficiente, 

Vom descrie în cele ce urmează, așa-zisele metode directe (datorate lui 
Gauss) adică metode prin care soluţia se obține într-un număr finit de pași 
(un număr finit de operații). 

„Apoi ne vom ocupa, pe scurt, de metodele. iterative (indirecte) în care 
soluția sistemului se obține printr-un proces de trecere la limită. 

Evident în metodele iterative nu putem efectua o infinitate de operaţii, 
deci ne oprim după calculul unui număr finit de aproximante (se face în acest 
caz o eroare de trunchiere). 

În metodele directe, dacă, s-ar lucra ideal (deci fără erori de rotunjire), 
“lucru imposibil dealtfel, s-ar obţine soluția exactă, 

| În timp ce în metodele directe erorile de rotunjire devin foarte periculoase, 
conducînd la. soluţii fără sens dacă sistemul este rău condiționat, în metodele 
iterative, în ciuda erorii de trunchiere erorile de rotunjire nu se acumulează. 


$ 2. Metoda eliminării a lui Gauss 


Pentru ilustrarea acestei metode directe vom utiliza un sistem de trei 
ecuații liniare algebrice cu trei necunoscute, despre care presupunem că admite 
soluţie unică: 

Aita AaYat dis = bi 
Ca Xi f- GaoVa + Aaaa = ba (5) 


(as Yi k aoa -F lasa = ba. 


Cel puţin unul din coeficienții aa, dz, da este diferit de zero, căci altfel 
sistemul (5) n-ar mai [i cu trei necunoscute. Dacă au, este egal cu 0 rearan- 
Jăm ecuaţiile astfel încît coeficientul lui 4, din prima ecuaţie să nu fie zero. 

Soluţia sistemului rămîne, evident, neschimbată dacă schimbăm ecuațiile 
între ele (De ce?). 

ua ară Clay 
Definim un mMultiplicator ma = — 
Au 


Înmulțim prima ecuaţie din sistem cu ma şi o scădem din ecuaţia a doua 
Ei : fi i ee Viei .. . E fs 
(atenție acum „prima și a „doua“ se referă la ecuaţiile deja rearanjate dacă 
a fost cazul). : 
După efectuarea operaţiei de scădere, ecuaţia a doua din sistem va fi: 
j: i | 


(da — Mada) A (22 — Pada3) a (23 — Mod43)X3 — 
— Do mah (6) 
Dar: 


a 
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du = 0 
Au 


Aa — Molly — Hal 


deci x, a fost eliminat din ecuația a doua [acum se observă de ce a fost 


da 
ales ma.= — |: 


dg2 = Aa2 — PMada2 
(3 = (as — Mada 
bi — ba — ?Mabi, 


astfel (6) devine: 


ta 
= 


Gina -F diaaa — bi. (7). 


Deci a doua ecuație din sistemul (5) va fi dată de relația (7). 
Soluția acestui nou sistem este aceeași cu soluția sistemului inițial. 


Similar, definim ma — Sti și înmulțim prima ecuație cu 743 și o scădem 
a 
din a treia. Eliminăm din nou coeficientul lui x, și obținem: 
Agatat gata = bg (8) 
unde: 
(ga — Aga — Misda 
33 = A33 — Palo 


bi — ba — mah. 


înlocuind în sistemul inițial a doua ecuație prin (7) și a treia ecuație 
prin (8) obținem sistemul: 


Aa FF aaa AF das Xa = bu 
Azaăta AF sa Xa — ba (9) 


Lă . ? Lă 
sa%a E 39 8 = bg. 


Sistemul obținut are aceeași soluție ca sistemul inițial și prezintă avan- 
tajul că x; nu mai apare în ultimele două ecuații. 

Să trecem acum la eliminarea lui x din ultimele două ecuaţii. 

Dacă asa = 0, permutăm ultimele două ecuații. Dacă și ai și ay sînt 
egali cu zero, am avea o infinitate de soluții sau nici una, ceea ce nu este 


"posibil deoarece am. admis că sistemul inițial are soluție unică. 


gasi Ass E 
Alegem multiplicatorul m = -%. 
ZE 


Dacă înmulţim a doua ecuaţie din sistemul (9) cu m; şi o scădem din a 
treia, obținem: 


(3 — m * âa) * Xa (aga — m = aia) - ma — be — ma: bi. 
Din nou dgz — 44 * Aa = 0, și notînd 
(33 — Aga — Ma * dia 
bg == ba — ma di 
obținem: 
53 ÎS, == (fi ! (10) 
Astfel putem înlocui sistemul (9) prin sistemul echivalent: 


Aaa Aaa%a A das %3 = bi 
Aaa Aaa %9 = ba (11) 
3 33 — by, 


Matricea asociată coeficienţilor necunoscutelor sistemului (11) este su- 
perior triunghiulară: 


L) 
du up 25| 
7 || 
U=||0 za All 
INI 
0 2 10) a || 


și de aceca sistemul (11) se numește superior triunghiular. 


Sistemul (11) se rezolvă foarte simplu, prin substituție, începînd cu 
ecuația a treia. 


Deci: 


bi — Aaa — (usXa 4 


du 


Amintim Că au 0, aia O iar ap este de asemenea, diferit de zero 
deoarece în caz contrar ar rezulta că matricea asociată coeficienţilor siste- 
mului (5) ar fi singulară (cu determinantul zero). 

În ecuaţia a doua şi a treia din (11) renotăm coeficienții CU dz, das, ba 
respectiv da, ba. Aceste ecuații se scriu: 


(ao%a At doza = ba 
Aaaa — ba. 


Deşi prin această renotare nu mai putem identifica pe ao2, da ba respectiv 
33, ba cu coeficienții ecuaţiei doi și trei din (5) renotarea este posibilă deoa- 
rece nu mai folosim niciodată acest sistem. Renotarea prezintă avantaj la 


scrierea programului de calcul pentru a nu introduce noi variabile ci folo-! 


sirea unor locaţii din program rămase libere. 

Operaţiile folosite în rezolvarea sistemului (5) asupra coeficienţilor și 
termenului liber nu modifică soluţia sistemului. 

Acum să trecem la descrierea generală a algoritmului lui Gauss. 

Să considerăm sistemul liniar de n ecuații cu n necunoscute, de asemenea 
presupunem că sistemul admile solutie umcă, deci că matricea coețicienților 
necunoscutelor este nesingulară: 


Aaa Aro fr eee Gap FF ee PF Da Xa = bu 
May ana ee da; Fl a di 
i RS pt MIO II ALE dn Rara rac aaa (12) 
Apr F Aia Da i a d) 
AX Aaa ! Aj) îi Ama bus 
Dacă notăm: 
As Aagece Mag == A] 
| 
|| Zap aa ie Go me aa] 
A | | 
| .| 
| manea LOT) se dau] 
| +] |2| 
| 
Î| ace | | 
„isi | Mm. =] ba 
| Sa | | : | 
i | 
| a | | ba | 
atunci sistemul (12) se poate scrie în forma matriceală,: 
i A ev SERIA (13) 


Presupunem că a 4 0. Eliminăm pe x, din ultimele (p—1) ecuații 
scăzînd din ecuația „1* prima ecuaţie multiplicată cu: 


[4 g 
Dau ot RS | Ei Za 2 
Fa 


deci avem p—l multiplicatori). 
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Li 
. 


Ultimele 4 — 1 ecuaţii nu vor mai conține necunoscuta x; vor deveni: 


o 
[i 


2 aţa, — 2) 


aţa) Aa A see lu aț?) Xp = b(2) 
unde: 


(2 
A) — App — du 


La) 


a ; 
69 — d; — mb 
D= di 


Acest sistem este un sistem de n—1 ecuaţii în necunoscutele %2, 3 dq... 
cpt 

Dacă at) z 0 similar putem elimina pe a din ultimele (p—2) ecuații 
în necunoscutele De uda 


n 
Fie: i 
Ma = sar) 
coeficienţii acestui sistem sînt al de: 
a) — af) — ma) 
V(8) == 002) = mmaabl?), 7 = 3, d cs n 


Elementele a, 47, a, ... ce apar în timpul procesului de eliminare a 
necunoscutelor poartă numele de elemente pivot. 

Dacă toate elementele pivot sînt diferite de zero, după n—1 pași obținem 
o singură ecuaţie: 


Ghid cs A) 


LO) 


Acum, colectind prima ecuație din fiecare pas obținem: 


AI aa a (9) Ata 2 „a te alt) bl) 
a) at al), = b2 


tel a Limite pl | li) 
3 


440) Xp = bt) 


unde am notat: a aș: bb = b, pentru coeficienții sistemului inițial pe 
care dorim să-l rezolvăni prin metoda Gauss. 


Am folosit în (14) o notație mai complicată pentru a arăta exact cum 
se transformă coeficienții și termenul liber în metoda lui Gauss. Remarcăm 
din nou, ca și în cazul sistemelor de trei ecuaţii, că la scrierea programului 
de calcul este mai simplu să revenim la vechea scriere a ecuaţiei. 


Sistemul (14) este un sistem triunghiular de forma: 


Gata A Aaa eee E Anania = bi 
Aaa ia iii om 1-1 iii Dap = ba 


Ap Xa- 1 + Anna = bai 


Am n =: by 


Dacă presupunem că a, 40, 1=— 1, 2, „n atunci necunoscutele x, 
„ %, se determină astfel: 


1 
ia 
DEE : : 
x, = (din ultima ecuație) 
dan 
e) Ca rusa . Ş 3 
e e tie nn n (din penultima ecuație) 
A, a-l (16) 
g = b, lin Xa Aa „n-1%n-a — Aaa 
sil 


Relaţiile (16) se pot serie compact astfel: 
ij => > did 
js Ea e e la my Îl (17) 


i 
a 


Matricea coeficienţilor sistemului (15) se wuwmește superior triungulară, 
În cazul în care am avea 0 matrice inferior triunghiulară, deci cu zerouri 
deasupra diagonalei principale, rezolvarea acestui sistem este aceeași numai 
că se vor scoate necunoscutele în ordinea Me Manea Am 

Dacă studiem formula (17) observăm că pentru rezolvarea unui sistem. 
triunghiular se fac „n“ împărțiri şi 

EiĂ, Pa l lee deo 
NR (4 — 1) zi: n(u — 1) adunări și înmulţiri. Ş 


| 
Dacă „n“ este „mare“ atunci n >" și 
1 | 


nl -— 1) a soti 


yu A = .. SE. 
Evaluăm în calcul numărul de operaţii pentru a ne da seama cît timp 
calculator necesită o anume problemă și deci cît costă rezolvarea ei. 


'$ 3. Formularea matriceală a algoritmului lui. Gauss 


Urmărind cu atenţie algoritmul lui Gauss se observă că operaţiile care 
se fac asupra coeficienţilor a; se fac și asupra componentelor vectorului 
format cu termenii liberi. A 
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Formăm matricea A& (extinsa matricii A) din matricea A adăugîndu-i 
o coloană formată din vectorul 2. 


Deci: 
e da du A b, | 
| Elea Do 08 07 deea Il (18) 
= [si Vlupiaae Che oo ireal e Ma 


Operaţiile din algoritmul Gauss ce se efectuează cu liniile din matricea 
A& sînt: 
— înmulțirea unei linii cu o constantă, 
— scăderea, unei linii din alta. și înlocuirea celei de a doua cu aul 
D + scăderii, 
— rearanjarea liniilor. 
Dacă notăm cu a; pe d, atunci matricea A& cu „n“ linii şi „n +1 
coloane este: 


Gai Aaa din Aumi]] 
Li] 
day Maas A2p e amr (19) 


Aaa ma am amr 


Wa aa e Aug LATA im 
0 Aaa ve Aa m-a dan Manu 

7 a Is | eo ada ee O Oi el e pa Asia (20) 
0 0 a i Anti An-Lmti 


0 Oz 0 dan Ami 


Se observă procedeil de renotare a elementelor matricii A& (vezi orga- 
nigrama- și programul) căci elementele a, din (20) se calculează după cum 

| s-a arătat în paragraful precedent. 
| În procesul de elimiriare se observă că la eliminarea necunoscutci x, ele- 
| mentele matricei din coloana £ și liniile & +1, ..., nl = 1, 2, sm, n — 1) devin 


| nule. 
Ş 4 ) luă do 
Înmulțim acum ultima linie din matricea (20) cu — și obţinem: 
: dan 
Ami 
( (0) as O al zi), (21) 
dun 
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PEPE 


Renotînd ultimul element al vectorului (21) cu aa obținem: 


(00040. 28 Da Iza a (22) 


Ultima. componentă a veetoruliză (22) este componenta +, a vectorului 
soluție. 


Înmulţim pe (22) cu au și scădem din linia „n— 1“ și astfel obținem: 


(0, 0, 0, neo al p—1 0, Ami — ao ua +1): (23 


Înmulțind vectorul -(23) cu obținem pe ultima poziție (deci 


a 
nun 
în coloana n+ | a matricii A&) valoarea necunoscutei x, a. 
Continuînd procedeul obținem pe rînd necunoscutele 3,25 Xp 35 ns dtz, Xa 
și matricea A* devine: : 


| aur 
000.10 Ea eat 
Co al 
At SU (24) 


N 


Q N e, | 
| 


Cami | 


Observaţie. — În matricea (24) elementele de pe ultima coloană sînt 


exprimate de relațiile (17) cu b, înlocuit cu a; psa: 


$ 4. Organigrama pentru metoda lui Gauss și programul de calcul 
în FORTRAN 


În paragrafele 2 și 3 am studiat metoda eliminării a lui Gauss și s-a arătat 
că pentru fiecare pas „7“ în procesul de eliminare este necesar să avem a, z 0. 
Dacă a, — 0 se caută în coloana 7 un element diferit de zero, lic acesta a, 
și se schimbă linia / cu linia 5. Dacă se foloseşte renotarea înseamnă că averi 
a; 40. 

În cursul acestui capitol ($ 6) se va arăta.că pentru a obţine o soluție 
a sistemului algebric liniar, cît mai apropiată de soluţia exactă, nu e necesar 
numai ca 4, £Oci trebuie să avem: 


a; = max la] 
Shu 
(elementul pivot să fie cel”mai mare, în modul, în coloana 4). 
Algoritmul lui Gauss de rezolvare a sistemelor liniare algebrice descris 
în paragraful 2 poate fi urmărit cu ușurință în organigrama din figura 11.2. 


Urmează programul de calcul corespunzător, în limbaj FORTRAN în 
care sînt demarcate prin cartele de comentariu blocurile de calcul din 
figura II. 2. 
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Rezolvă sistemul cu matricea asociată | 
da coeficienţilor necunoscutelor superior triunghiu- 
lară 


nu 


. 
- = a tipărește soluţia 
Aflarea elementului maxim | în modul) al 


| coloanei „i. Fie „ k” linia acestui element. 


Este elementul maxim pi felu fete Matricea este 
in modul) egal cu zero ? singulară 
> E STOP 


Schimbă linia „ka elementului maxim cu 
lima „i Aceeaşi operație cu termenul liber 
(schimbă by cu b;) 


Împarte elemetele liniei „i ” şi componenta bi 
cu elementul pivot aii 


Scade linia „i ” înmulțită cu ajj din „linia „j” 
1V Scade din bj, pe bi-aij 
| Solia do pomcugyl) 


Fig. I1.2. Organigrama pentru metoda eliminării a lui Gauss 


PROGRAMUL FORITRAN 


peniru metoda ehmunării a lui Gauss 


SUBROUTINE GAUSS (A,B,N,EPS) 
DIMENSION A(N,N), B(N), 


AICI ÎNCEPE PROCESUL DE ELIMINARE 


CL eio 


DO 1001 =1,N 
AMAX — A(LI) 
IMAX —1 


AFLAREA ELEMENTULUI MAXIM DIN COLOANA | 


(DONE) 


DO 101] =1,N 
IF(ABS(AMAX) — ABS(A(J,D))).102, 101, 101 


102 AMAK = A(J,D 
IMAX = ] 
101 CONTINUE 


(e) 
C VERIFICĂ DACĂ DETERMINANTUL ESTE DIFERIT DE ZERO 
(e; 


IF(ABS(AMAX) — EPS) 103, 104, 104 
WRITE (108,3) 

FORMAT (1X, "SISTEM SINGULAR”) 
STOP 


10 


e) 
3 


C 
C PERMUTAREA DE LINII, DACĂ ESTE NECESARĂ 
(6 


104 DO 105K=I1,N 
T = A(IMAX, N) 
AUMAX, K) = ACR) 
105 A(LIK) = TJAMAX 
T = B(IMAX) 4 Î 
BUMAX) = B() 
B(1) — T/AMAX 


c 
C URMEAZĂ ELIMINAREA NECUNOSCUTEI S(1) A 
Cc 

TF(I—N) 106, 100, 100 


AG ip de-ai 
DO 107 ] — JJ,N 


XM = A(J,D 4 
DO 108 E =1,N 
108 ACR) = A(R) — A(LR)-XM 
107 B(J) = BD) = BO)-XM 
100 CONTINUE ; 
ă | 
C AFLAREA SOLUȚIEI PRIN SUBSTITUŢIE 
c 


DO 109] =2, N 
E=N-—J+a 


ÎI Dee ȘI 
DO 101 =1,K 
110 BUD = BO) — AQ, KE) - B(KR) 
109 CONTINUE 
E RETURN 


END 


Exemplu numeric 
Fie sistemul: 


What is=l 
Xa 10001 x3-+ 23 = 2 | 


pe p> Poeta a dea = I| 


pe care dorim să-l vezolvăm prin metoda Gauss. 
La pasul 1 în procesul de eliminare avem: - | 
— elementul maxim în modul în coloana 1 este |, deci elementul pivot 

SE Il e Îl 
- întrucît elementul pivot este egal cu -unitatea, prima ecuație a siste- 

mului este deja normată (împărțită la elementul pivot), 
— eliminarea necunoscutei x, din ecuația a doua (se scade prima ecuație 

din a doua, deoarece a = 1) ne conduce la: 


0,0001 4 +aăs=l, 
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— eliminarea necunoscutei +, din ecuaţia a treia ne conduce la: 


Xa ăg = 0. 
Sistemul de ecuaţii devine: 
Lai Xa 3 = | 
0 + 0,0001 += 
| 0 + cr leat), 


La pasul 2 al eliminării avem: ca 
Aaa = 0,0001; az2 = |, şi elementul pivot este aa» (este element maxim 
în coloana 2), 4 ! 
— se schimbă linia 2 cu linia 3 şi sistemul devine: 
Xa a = | 
Xa Xa = 0 
0,000lx3-+ a = | 


— eliminăm necunoscuta x2 din ecuaţia a treia și avem: 


Mat iata l 
Xp Xa =0 
0,9999x = 1, 


La pasul 3 urmează rezolvarea sistemului cu matrice superior triunghiu- 
lară (se operează cu trei zecimale) 


d Xa— sl (—1) 1 Bl 
Același sistem, rezolvat cu calculațorul, a furnizat soluția: 
3 = 1000000; wa — —1,000099; xs — 1,000099 


(calculul s-a efectuat cu 7 zecimale). 


$ 5. Strategia pivoțării 


În descrierea algoritmului lui Gauss noi am presupus (vezi $2) că elemen- 
tele aţ% (elementele pivot) sînt diferite de zero, altfe] sisțemul este singular. 
Să considerăm următorul sistem de ecuatii: 
Wa Xa = 
Xa Xa 25 = 2 (25) 
Wa dak dis = | 


Sistemul admite soluție unică (determinantul matricei coeficienţilor este 
nenul) xi =, = —l, xa=l. : 
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Aplicînd procedeul lui Gauss, după primul pas primim: ce 
ina =ll 
Xa Ay = 0. 
Deci, observăm că a = 0 (coeficientul lui x> în prima ecuație (26) — 
vezi $ 2). 


Pentru a putea ae acest neajuns permutăm cele două ecuații și 
scriem sistemul (26) c 


(26) 


Xa-ţ- 33 = 0 


Bin 22) 


şi acum obținem un sistem deja triunghiular. care rezolvat ne conduce la 
soluția căutată. 

Să considerăm acum cazul Bea Deci lapasul £ am obținut ap = 
Atunci măcar unul din elementele aţ a pa Va i „aţi este nenul, A caca 
cea A ar fi singulară. 

Să notăm acel element cu a4.z 0. Vom a atunci la schimbarea 
între ele a liniilor k și 7 apoi vom continua eliminarea, 

Deci orice sistem nesingular poăte fi! redus la o formă triunghiulară cu 
ajutorul metodei lui Gauss combinată, cu schimbări de linii între ele. 

Să revenim. acum la sistemul (25) în care coeficientul lui x2 din ecuația 
a doua (425) îl înlocuim cu 1,0001. Rezolvînd din nou prin eliminare obținem: 


Xub at ăp=l : i 
0,0001 %34- x = | (28) 
—9 999x3 — — 10 000. 


Acum, rezolvînd sistemul (28) ca un sistem triunghiular și utilizînd trei 
zecimale obţinem: 


aa 2 0imsza =— Oiaral 201,000 (29) 


Am scris 1,000 pentru a scoate în evidenţă că lucrăm cu trei zecimale, 
Soluţia exactă cu rotunjire la patru zecimale este: : 


ata = 1,000, aa = 1,0001, xp = 10001. 


Acest neajuns provine din faptul că aţi) a devenit după primul pas al 
eliminării 0,0001 (deci zero cînd lucrăm în “trei: zecimale), așadar:un element 
pivot nu- are voie să fie nici apropiat de zero. 


Noi am rezolvat în $ 4 sistemul: 


aka al 
Xu 1,000 la 23 2 2 
sut Dă 23 = | 


prin schimbarea liniilor 2 și 3 şi am obținut soluția: 
Xy > — a = Ba — 1,000 


corectă. în trei zecimale (am utilizat trei zecimale în calculul soluției). 
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Exemplul piezentat nic arată că erotile de: rotunjire pot avea efect ca 
tastrofal în metoda Gauss, de aceea foarte des la pi k urmăm strategia 
(numită strategie de pivotare): 


— căutăm cel mai mic număr pr : Dentru care: 
Jal?| = max la] ksisu 


Și FD păzi) liniile k și 7. 


Cu alte cuvinte, pivotul la fiecare pas este cel mai mare "elertent în modul 
în: coloană '-considerată. : 

Deci, începem prin. a căuta pivotul în prima. coloană. Linia întîi și linia 
căre conţine „elementul, pivot (cel mai mare element. în modul „din coloana 
întîi) se schimbă între ele. 

Procesul se repetă cu a doua coloană și așa inai departe. 

Cînd linia ce conţine pivotul este înmulțită cu un coeficient (înainte de a 
o scădea din altă linie) ea este de fapt înmulțită cu un număr subunitar 
(vezi definiția multiplicatorilor 4, $ 2) și astfel erorile de. rotunjire în liniile 
pivot sînt înmulțite cu un număr mai'mic ca |. 


Observaţie. Procesul eliminării a lui: Gauss are loc fără schim- 
bări de linii între ele dacă matricea coeficienţilor necunoscutelor 
A > alia 


=> ta 


are proprietatea că dă diagonal” dominantă, adică dacă: 


deci elementele de pe'diagonăla principală sînt mai mari (sau ega- 
le) în rolul „decâț, suma, medulelor. elementelor. din linia respectivă. 


Deexemplu: . îi mai 
| pi 220 aj 
2 = 


este o matrice diagonal dominantă de ordinul trei. 


Exerciţiu. Este matricea de:'ordinil 1x37 'diagonal. dominantă ? +, 


MI 32 ii e to 10 o jpâeu-0 
ri 2iitii foto 0 - :0 

DRE 20 1 Pda! +6 

De ni lepra o ae 90 ari 0 
i pi Rac A Aa pp det at ue 29 a Ad) 

OR OR 0 NORA = 


1 


0: 0 :0 0 tr0—l 2 


(Astfel de matrice se numesc tridiagonale și apar foarte frecvent în practică, 
Se numesc matrice rare, avînd multe elemente egale cu zero.) 
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| $ 6. Sisteme rău condiţionate şi erorile de rotunjire Alla = max (10%. (0,1441 + 0,8648), 108 (0,2161 + 1,2969); — | 


ÎN 1 Mala | 
|N Dacă considerăm sistemul Ax = b unde: . : 2410 II | 
| | id at, auz ||. mă x A: d, lar: | 
|| (oi aa j Xa i ba Ala = 2,1617. | 
și Z soluţia sa calculată aproximativ, vom numi vector rezidual r, vectorul Deci 
|] AR DI ae Î conda(4) = 2,1617 - 1,5130 - 108 2 3,3 - 10%. (32) 
| E da ajăudi Pest6, solii jiag ea ct „cleci gi aaa iat ei a Relaţia (32) ne arată că sistemul considerat este foarte rău condiționat, | 
| aștepta că dacă cele două componente ale vectorului rezidual ar fi mici, stea i (m Dau alu c ie old vaaticieţie utira bcu opace: Satan ara e ul 
atunci Z să fie o soluţie suficient de bună pentru sistemul considerat. Asia Ie 
Lucrurile nu stau nici pe departe astfel: Observaţie importantă. Chiar dacă elementele matricei A și 
Iată următorul exemplu foarte sugestiv: a vectorului b sînt date exact, totuși din reprezentarea în calcu- 
| d Fie: lator (virgulă mobilă) rezultă o matrice rotunjită A cu elemen- | 
| a dn | 1,2969 0,8648 tele â,, unde: ij | 
| 0,2161 chiu | | Zeta 
II , || 0.86% . la] | 
| | "| o;14240 | (4) u reprezentînd unitatea rotunjită în reprezentarea numerelor în virgulă | 
|| A mobilă pe calculator. În cele mai multe mașini de calcul „4“ se află între | 
| Presupunem că: i | 104% și 10-6. | 
N | i | 
| p = | îi] , | $ 7. Metode iterative pentru rezolvarea sistemelor liniare. | 
ij [| 0,870] Metoda Gauss-Seidel | 
| . : E: 
| Atunci: | Să considerăm următorul sistem de două ecuaţii cu două necunoscute | 
Wu 7 0,8642 1,2969 0,8648 0,9911 —10% . : SE | 
Vl pi || [RE | — i -] = a SĂ 9 - (34) | 
72|] [| 0,1440 0,2161 0,1341] || —0,4870 108 PA Mai e 
| Deci componentele lui 7 sînt destui de mici, de ordinul 108. Totuși Z cu “Soluţia sistemului este x=1, y=1. | 
componentele daţe este foarte departe de soluţia exactă, care este: „Să | 
2 | 
| o = | |. | 
—3] | 
| 
| 


Sistemul considerat este foarte rău condiționat. După eliminarea lui x, 
obținem: i | 


a, = bf), unde: | 


0,2161 


| aţz) — 0,1441 — -0,8648 — 0,1441 — 0,1440999923 2 103. 


IN] . A 
| 
[| Deci o mică schimbare în coeficientul 0,1441 provoacă o mare schimbare 
în a, deci în x. i 
Dacă inversăm matricea A, obținem: 
4-1 pa] 01441 —0,8648| 
Dl 0261 1,2969 | 


Scoatem din prima ecuaţie pe x şi din a doua pe y. Obţinem: 


pe, (= 3) 


'8, 
(35) 
1 fi 
DE pa (o). 
| al da 7) 
Vom considera următorul proces iterativ: 
| 
x = —(4 — pd 
papi ai) 
(36) 
| 
ACE) ci a (7 20) 
elita) 


plecînd de la o aproximaţie inițială (de regulă se pleacă de la x == 0, y=0). 
„ Deci fie w0 = 0 prima aproximaţie, atunci din prima ecuaţie a lui (36) 
obținem: 
j | 4 i ei AR E 
203) dea zi 4= E (vezi graficul din figura 11.3), 


„Apoi, din a doua ecuație: 


4 10 
pb) | 2 i = (figura II.3) 


4 A 


Continuînd procesul de calcul din (36) scoatem: ii 


e [a oaze, 


3 9 pi ? 


N [ . că 


3 21] 3 


Observăm că deși am făcut numai două iterații ne-am apropiat destul 
de mult de soluția sistemului (se vede de asemenea și pe figura 3). 
Dacă scriem matricea asociată coeficienţilor necunoscutelor sistemu- 


lui (34): 
8 1 | 
1 —3 
observăm următoarele lucruri: 
a) elementele de pe diagonala principală sînt nenule și mai mari în modul 
decît celelalte elemente. 


d) aul > lana] = 3 lan] = 3 lazi], 


deci elementele de pe diagonală domină cu raportul trei celelalte elemente. 


Geometric, observăm că panta primei drepte (corespunzătoare primei 
ecuații) este —3 (în modul mai mare.ca 1) iar a celei de a doua este mai 
mică decît 1. 
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Acum, să vedem ce se întîmplă dacă permutăm ecuaţiile între ele, Deci: 


pa > ap) (37) 
Va 4— 3%. 
Pornind cu acecași aproximație inițială (4 = 0, y = 0), obținen» 
Iterația x y ; 

1) 0 (0) 

1 —2 10 (38) 

Zi 28 —380 

3 — 242 730 


Fig. IL. 4 


Ce s-a întîmplat ? 

A . 

Panta: primei. drepte este mai mică decît 1 iar a celei de a doua este mai 
mare decît 1 (în modul); : 

Urmează că pentru stabilirea convergenței va tvebui să studiem matricea 
coeficientilor necunosculelor sistemului. | 

Considerăm pentru aceasta. forma. generală a unui sistem de două ecuaţii 

cu două necunoscute: 


5 r 


At day = bi a 1 (39) 
Aaa daay = ba. : 


'39 


Procesul iterativ este definit prin: 


A 1 

EA => că (ip ay) 
du (40) 

în). | b (le) 

5! opel a ). 
22 

Notăm: 
Ab = pp — x 


Ay = y— mb, 
Din prima ecuaţie (39) și prima ecuaţie din (40) deducem: 


dia 


At) SS Awt-b 


du 
ȘI similar, din celelalte două ecuaţii: 


Ay = — da A ah). 
Aaa 


E] 


Din ultimele două relații deducem: 


: Ma" a 
Ay) — a Cai Ayk-v, (41) 
da * dap 
Analog, 
Ay = Aia? da 4 pa, 
Dat dur * Aa 
eci: 
. 2 
Ai = pa) A p=2), 
Aa * Aaa 
Continuînd raționamentul, obținem: 
Ş di-a 3 z Lă 
Ay = [E 22) Agia) =, — [Zi Za At 42 
Ai * Aaa du" A ai ai, 
aL? Ma 
La tel se obține: 
. Li 
R) __ [| Bia * day 
Ati) — (2 =] Azto, = (43) 
di * Baa 


Relaţiile (42) și (43) ne spun că șirul i k, ă 

i i „Spun că șirul de iterate (x) converge către x 
(prima componentă a soluției sistemului (39)) și cd, aici Ste y (a 
doua componentă a soluției sistemului (39)) dacă: 


Aaa 
—a 21 E ll (44) 


Zar * Aaa 
Condiţia (44) este îndeplinită” dacă, de exemplu: 
[iza > | aaa] 

aa] > lazi] 


40 


sau: - 
(ist > lau] (46) 

laz2] > azil. 
Cu alte cuvinte am ajuns la concluzia că diagonala matricei asociate 


sistemului domină celelalte elemente. i a 
Vom descrie în continuare metoda Gauss-Seidel pentru sisteme liniare 


de tre: ecuaţii cu trei necunoscute. 
Fie sistemul: 


isa Aaa Xa F das Xa = bu 
dot aa Xa t da 3 = ba 
daf Aga Xa F dag 3 = ba. 
Urmează că sistemul (47) se 


(47) 


Presupunem că du 40, da 40, as 70. 
poate scrie sub forma: 


ji = (d. Aaa 13%) 
un 
: 1 
iz = (ba AA (2p%3) (48) 
d2a 
| 
Xa = — (ba — Aa — 322), 
433 


Vom nota prima aproximaţie cu x(9, 40, a40). 
Calculăm noile componente după cum urmează: 


1 
34 = — (ba — az) — d13309)) 
Au 
1 
aţi) (Ba — az Aag40)) (49) 
Aaa 
1 
i = (bs — dsu%( — asa), 
433 


Astfel s-a încheiat prima iterație. Acum punctul 4, aţi), „x devine 
aproximaţie inițială pentru iteraţia a doua (ia locul, în memoria calculato- 
rului, punctului 440), 40), 240). 

În general algoritmul este dat de: 


i 


a — — (bn — aaa D — apa) 
du i 
(5) 1 : -(R) „(1 50 
A => — (Ba — daf) — app 1) (50) 
A2a2 
(E) 1 (i) pla) 
a (Bas — asa) 32%), 
433 
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Putem trece acum ușor la algoritmul Gauss-Seidel pentru un sistem 


de n ecuații cu n necunoscute, 4 ' 
Dacă matricea coeficienţilor este A = |a,,|i<î 
91 
bi | 
iar vectorul termenilor liberi este: b=— | : | atunci: 
| 
i n j | 
PD ra O ee a EU 

afet = Fi e (60) 


du 


da 


iesionaast la 2 Se 4 | 


Ca primă aproximaţie se consideră, de obicei, un vector n-dimensional” 


cu toate componentele nule. E 
SUM=SUM57 ALL) XI) 


Calculul se opreşte cînd: | ia 
oase acţ) made), <a les (61) | 
Let 4 

Altfel spus, cînd componentele vectorului soluţie la două iterații succe- | 


sive devin foarte apropiate înseamnă că ne-am apropiat de soluţie. | 
Din discuţiile anterioare se poate deduce că metoda descrisă converge 7 
cînd: | 


lau] > lalea m lau (62) 


SUM = SUM+S 7 A(1,J).X01J) 
J=1 j 


a : . X1(1) = (B(L) -SUM)/AL(I,L) 
pentru orice 1 = Î, 2, 3,...., n şi pentru cel puțin un 4 avem inegalitate strictă. 
Deci: 

ER =maxlER,|XI(1)- X01I)|) 


lau] > leat aa st leul. (6) | | 
De exemplu pentru matricea nxn 
De al : 
—1 2 —l 9] 
—l 2 il | £ | 


—1 2 —l » 
4 = (64) 


Q i 2 | 
metoda Gauss-Seidel converge. 

X0(K) = XI (K) 
Exerciţiu. — E xplicaţi de ce! Set as 
Dăm în continuare, organigrama (fig. II. 5) și programul de calcul pentru 


metoda Gauss-Seidel. Urmează de asemenea un exemplu de calcul. si, 15. Osea taia adie iza PO taie Craig Se] del 
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Programul FORTRAN 
pentru metoda Gauss-Seidel 
SUBROUTINE SEIDEL (A, B, X0, XI,N,EPS) 
DIMENSION A(N, N), B(N), XO(N), XI(X) 


GI 
C  VECTORUL X0 ESTE APROXIMAȚIA INIȚIALĂ, Î ) SE E 
C SOLUŢIA SISTEMULUI A IA zu Nat 10997 SI oa 
(e 
8 ER =0 
îi il 
4 SUM = 0 
E aa i! | 
îi) il 
1F(10)2,2,3 
a DO 10 ] = 1,10 | 
SUM = SUM + A(I, ].X1 
+ TR(UL—N)2,2,12 a) 
2 130) pi ET 
i SUM = SUM + A(I,]).X0(J) | 
12 X1() = (B(D) — SUM)/A(I) | 


ERI = ABS(X1(1) — X0()) 
ER = AMAXI(ER, ERI) | 


I=I1+1 4 
IF(1—N)4,4,5 , | 
5 IF(ER — EPS)7,7,6 
6 DO 20K = 1,N 
20 X0(I) = XI1(1S) 
GO TO 8 
7 RETURN 1] 
END 


Exemplu numeric. 
Să. se rezolve, prin metoda Gauss-Seidel sistemul : 


A-x=—b, unde: 
| 4 =1-0 1 | 
—1 4.0 — | 
A = 5 l b = 2 . 
—] 0 4 —L|P 1) | 
0 -—1 —l 4 | 
Dăm mai jos tabelul cu primele 10 iterații în rezolvarea acestui sistem: $ 
Nr. iterației Zi EP a A ă 
i] 0) 1) 0 1) 
1 0,25 0,5625 0,0625 0,40625 
2 0,40625 0,703125 0,203125 0,476563 
3 0,476563 0,738281 0,238281 0,494141 
4 0,494141 0,747070 0,247070 0,498535 
| 5 0,498535 0,749267 0,249267 0,499634 | 
6 0,499634 0,749817 0,2498 17 0,499908 | 
7 0,499908 0,749954 0,249954 0,499971 N | 
8 0,499971 0,749989 0,249989 0,499994 | 
9 0,499994 0,749997 0,249997 0,499998 | 
10 0,499098 0,749999 0,249999 0,499999 | 


Soluţia exactă este: (0,5; 0,75; 0,25; 0,5) şi în cele 10 iterații s-a atins | 
precizia de 1,5 : 10% : 
Exerciţii 

1. Rezolvarea sistemului de ecuații: 


(aer Aaaăa = bi 


Gas t Aaaa = ba 


tul de intersecţie al celor două drepte date de ecuaţiile 


vevine la a găsi punc 
le două drepte verifică inegali- 


sistemului. Să se arate că unghiul făcut de ce 
tatea: 
i / i N 
|cts a| < să cond, (4) (se presupune matricea A  nesingulară). 
2. Arătaţi că printre toate matricele 2x 2, nesingulare, ale căror elemente 
sînt întregi nenegativi şi nedepășind numărul 100, matricea 


| 100 99 
| 


99 98 


are numărul de condiţionare euclidian maximal. 
3. Să se calculeze cond„(4) pentru matricea asociată sistemului: 


Sea Slide) 
6x — 3,979 — 2,03. 
Cum se schimbă soluția acestiii sistem dacă membrul drept devine: 


pre i) ] (deci s-a schimbat puțin membrul drept). 


4. Este nevoie la rezolvarea sistemului: Ax = b unde: 


| ei il o | 
faza ta o la milă 
dica poe caen el Atu] 
eee 4 1 || 


de schimbarea liniilor între ele? Explicaţi răspunsul. 
5. Rezolvaţi, pe calculator, prin metoda lui Gauss sistemul: 
9x, — 2x34 Xs = 50 
i e Sita — 3%a = 18 
— 2 + dia 13 = 19 


utilizaţi în calcul 6 poziţii zecimale. 
6. Scrieţi un program FORTRAN pentru metoda Gauss-Seidel și rezol- 


vaţi prin această metodă sistemul: Ax — b, unde: 


4 —i —l 0], | 
Se both dies 2 
Ala 0. A ZU b=| 9 


a 


0 eilpez ul | 

a) În câţi paşi obţineţi o bună aproximaţie? (lucraţi cu 5 cifre zecimale). 

;) Comparaţi rezultatele c& rezultatele obținute prin metoda Gauss. 

7. Să se dea o margine inferioară pentru numărul de condiţionare 
cond, (4) unde: 


| 1 —l 1 | 
Azul sal e el, ez 0 
] 1 e s| 


i 8. Evaluaţi perturbarea eventuală a soluţiilor sistemului: 
x—2y = —1l 
|—2x + 40ly = 2 ! 


dacă, schimbăm componentele membrului drept cu 0,01. Găsiţi soluţia siste- 
mului cu aceeași matrice a coeficienţilor și cu membrul drept dat de: 


5 = | pd | 
201 | 
9. Arătaţi că schimbările de linii :sau de coloane într-o lila A nu 
schimbă numerele de condiționare cond, (4); cond, (4). 
Exerciţii facultahve (cap. II ) € 
Se dă matricea (de dimensiune n x) 


d Sa SOMA mu Hi 
„peeestitoze 000 a ea za 0| 
E ea e = 
pi cedeaza i să] 


a) arătați că determinantul lui A este n +1; 
b) rezolvaţi ecuația de gradul n: 


3 


2—A —l 0) e Asa eta 0 

—1 De A Dog alle pe ea ati 0 

E oi Cora dez Atac eo pa AU [| 200) 

0) sola oare oaia e iei 2 

0. PI SE 03 2 pa 2 A 
utilizînd substituția A == 2 cos 8. i 


__£) Notînd cea mai mare rădăcină a ecuației cu A și cea mai mică. cu 
A) calculați: Sia ; 3 20 
„E Al) 

Îbhnay Bă 


1% N) 


d) Arătaţi că inversa matricei A este dată de: 


(fc ella) e u) v 
Sia stie i 
n+lg i 
Au 2 
o 0 e ie pa 
' | i 


e) Aplicați algoritmul lui Gauss pentru a inversa matricea A, de dimensiuni 
10xX 10, utilizînd în membrul drept matricea unițate. 
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Capitolul III 
REZOLVAREA ECUAȚIILOR NELINIARE 


$ 1. Introducere. Metode iterative 


Dacă f(x) este o funcţie reală de variabilă reală, atunci soluţiile ecuaţiei 


f(x) = 0 în general nu pot fi exprimate sub o formă convenabilă (0 formulă 


de exemplu) și în consecință, pentru obţinerea lor, se adoptă metode itera- 
AA aa aia 

„Pentru a ne face o idee despre ceea ce înseamnă, metodă iterativă, în 
acest cadru, vom considera ecuația: i 


ar) În E pati 


unde F(x) este o funcţie (de regulă derivabilă de mai multe ori) ale cărei 
valori le putem calcula pentru x aparținînd unui interval. 
Se pleacă de la un punct x, aproximaţie inițială, și caiculăm șirul: 


= Faro) sasa 0000 ))5 oda zu le (072 em (2) 


Fiecare etapă din calcul: 


ru = F(x) (3), 


se numește iteraţie. 
Dacă, şirul (x, converge către o valoare limită a atunci 'din continuitatea 
funcției [(x) rezultă că: : A 
lim F(x) = Fo) 


deci, trecînd la limită în (3) obținem: 
a = F(o) Boa 8 ți ide) (4) 


adică « satisface ecuaţia x = F(x). A i 
O interpretare geometrică a acestui calcul iterativ este dată în figura 111.1. 
Se vede că, plecînd din 4, obținem F(xo). Apoi ducem o paralelă la axa 
Ox prin punctul (40, F(40)) pînă intersectează prima bisectoare y = %, cobo- 
rîm din acel punct perpendiculara pe OX și obţinem pe xu(4 = F(40)), apoi 
Xa = F(x.) și procedeul continuă, 


47, 


Fig. TUI. 


Dacă urmărim graficul funcției din figura III. observăm că în apropierea 
soluției pantele tangentelor la graficul funcției y — F(x) devin foarte mici 
(pozitive). 

Acest fapi ne conduce la a considera comportarea primei derivate în 
vecinătatea soluţiei. 


Utilizînd teorema de medie a lui Lagrange obţinem: 


Ata Tm Fl) — Paza) Eta) (5) 


unde £, se află între x, și x. 

În (5) am utilizat și relația de iterare x = F(x). 

Urmează. că, dacă |F'(x)| < 1 pentru toţi = într-o vecinătate a rădăcinii 
care conține pe xo şi x, atunci avem convergența asigurată, 

Acum prezentăm riguros demonstrația următoare: 

Teoremă. Dacă ecuaţia + = F(x) are o rădăcină „a“ și F'(x) există în 
intervalul J = (x; | — a| < e) şi satisface condiția: 

|F'(2)| < m < 1 atunci pentru orice xp e J (deci orice aproximaţie îni- 
ţială în intervalul /) avem: 


2) x e Pe 0 2 e, 


5) Iti 2 E ca, 
n 


€) Singura soluţie a ecuaţiei în intervalul / este «. 

Demonstraţie. (a) Demonstrăm prin inducţie. Fie xy e J, presupunem că 
Xp-1 € J, atunci arătăm că şi x, e J 

Într-adevăr, conform teoremei lui Lagrange avem: 


da 4 Di (i) zi F(a) iz F'(0,) (X-a î a) 
unde t,e ]. 
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De unde: 


EA -. a] =F' (£,) (m-a a a) | Ss m a |caa u PT a| s mp 


deoarece 4, € ]. 
Urmează că x, e ] și (a) este demonstrată. 
(b) Avem: 


E = a.| Ss ma — al 


ontinuînd raționamentul obţinem: 
|, i i a | < m Fara) ui F(a) | > pie FE) (pe = a) | Ss m E 


m-a O lo 
(e J). 
De unde, recursiv se obţine: 


| Nae A [EA 
ez cal] ses re i cap ai 


Cum 0<m< 1 avem lim m —0. 


n-o 


Trecînd la limită în ambii termeni ai inegalității precedente obținem: 


lim |, — a| = 0, deci lim 4, = a, 


ue n 


(c) Presupunem că ecuaţia x = F(4) ar avea în J] şi rădăcina A za. 


Atunci: 


B= F(8) deci, 
a— B= Flo) —F(R)=F'(m)(a—bp); ne] 
de unde: 
ja — 8] < m-la—B|<la—8| 


deci o contradicție și implicația (c) este adevărată. 


$ 2. O metodă rapidă pentru extragerea rădăcinii pătrate 


Ecuația «2 = c poate fi scrisă și sub forma x — F(x) unde: 


(a) = (a )) c>0. (6) 


Soluţia ecuației x = F(3) este «=—yc (valoarea limită a algoritmului 
Xa = F(x) (fig. 111.2). 
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i Urmează, programul FORTRAN pentru calculul rădăcinii pătrate: 


1 FORMAT (GE 12.5) | 
READ (105 1) C,EPS 
XI = 1:5 | 

12 XI1 = (XI + C/XI)/2.0 | 
IF(ABS(XIL—XI) — EPS) 10, 0, 1 | 
Îi EI = XA | 
GO TO 2 : | | 
| 


y 10 WRITE (1082) C, XII 
: 2 FORMAT (8X, "RĂDĂCINA PĂTRATĂ A AN i E DAIILI | 
i "E 12.5, ESTE", E.12. 2.3) SA | 
STOP = e 3 LR | 
END || 


« « Dăm în continuare: tabelul cu primele trei iterații în calculul rădăcinii 
j ț pătrate din numărul .2: 
i 1 
Nr. iterației a Via = —— [au 4 N, | 
) i i 3 2 3 
Fig. III. 2 | pază | 
- Fig - 0, 1) : = 
1 1,5 1,416667 
: 3 4 se: Ar 2 1,416667 1,414216 
De data aceasta algoritmul se reduce la șirul inductiv: | 3 1414216 :  Vai4z14 | 
: 1 A Ă i (2 —.1,414213562...) | 
aere | Xa = |” (7) | 
2 Xp | 
Dăm în figura IIL.3 organigrama apoi, programul FORTRAN pentru tă $ 3. retea pari ÎmrruA j 


calculul rădăcinii pătrate a unui număr. 


De asemenea dăm tabelul cu iterațiile în calculul rădăcinii pătrate a lui 


2 plecînd de la xo = L,5. ÎN Calculul valorilor funcției F(x) este afectat de erori (de exemplu de erori 4 
( 
L 


de rotunjire). 


Vom face, în cele ce urmează, o evaluare a erorii după. un număr finit de 
iterații, 


Să notăm cu Zi, Za, Za cu Ba ee șirul valorilor calculate. De asemenea, | 
| 
notăm cu 8, eroarea din calculul lui PE (200) 9 


| 

| 

| Avem: 

îi Zu = F(Z) + Ba! 1 = 0, 4, 2, us n | (8) 

| Fie a rădăcina ecuației = F(x), deci: | 
- = F(a). (9) | | 

Tipăreşte 
soluţia E = | Scăzînd relația (9) din (8) obţinem: ! i! 
Xe] Î | 
Zau == F(z,) îi F(a) i ba (10) I| 


Aplicînd teorema lui Lagrange, obținem: 


E , i | 
; Ce e d 19 (7) (2, Fa a) al 3, (11) | 
| Fig. WII. 3. Organigrama pentru calculul rădăcinii pătrate. Ă ş unde să aparține intervalului determinat de 2 Și de &, 


Relaţia (11) se mai scrie: 
(1 E F'(E,)) (Za fiii a) a F'(3,) (2, a 7) île ap (12) 


Presupunem că |F'(6,)| «ml! și [8] s 8 atunci (12) ne dă: 


m i 
[Za — al [aa Zar 8 (13) 
| — m lL—m 
Primul termen al membrului drept estimează eroarea de trunchiere iar 
al doilea termen estimează eroarea de calcul. 


E depinde de eroarea 5, numai din ultima itera- 
— 7 
ţie şi deci 2, poate fi considerată o cantitate exactă cînd putem estima pe &,. 
Acest lucru ne conduce la următoarea concluzie: Nu este necesar să cal- 
-culăm cu mare precizie primele iterate deoarece rotunjirea în acește iterații 
nu are efect asupra rezultatului final. 
Se spune că metodele iterative se autocorectează. 
Studiind relaţia (13) se observă că pentru „n“ suficient de mare în dife- 


Eroarea de calcul 


rența | Zsa — a| rolul dominant îl are eroarea de calcul a Așa că, dacă 
— m 


după un „n“ suficient de mare vom continua iteraţiile, vom observa că dife- 


rența |Z,. — «| se comportă ca o mărime de ordinul —— . 
| — 


Exemplu. Pentru calculul lui (2 plecînd de la xp — 1,5, obținem: o 


2, = al. (13 | 2 ) = 1,415 (rotunjit la 3 zecimale), 


2 ap a 
z, i (+ | =) = 1,4142 (rotunjit la 4 zecimale), 
WM 


| PĂ ep î d 
D= — (+ + = = 1,41422 (rotunjit la 5 zecimale), 
Xa ; 


1 2 = a 
Z = — (+ E = | — 1,41422 (rotunjit la 5 zecimale), 
)) X3 
Pentru 2 aparținind intervalului generat de 2, și /2 avem: 
- R'(ţ, 1 
E a e 
1—1F(0)| 1—1F(4)| 
însă: 
“Fi 40)) = Rae. și observăm că 
Pra 


F'(2) > F'(%) pentru &, cuprins în intervalul (Za, J2]. 
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Deci: 


F'(z 1 
Hi su Za — 2] + : -3. 
1 —F'(2)| : 1—1F'(z)| 
Eifectuînd calculele avem: 
F'(24) = 0,00009 x 1» Iza 
1 


1 = F'(Za) 


|2a— 42] < 


= 1,0001, 


1 
e ea Dai escu aul Os 
|Za 3] 2 


3 s 1 10€ (precizia în calcul) 
şi cu aceasta: 
|24— VZ| <0;50005 - 10-10 1,0001 : 10”€. 
Se observă că după 4 iterații precizia depinde numai de eroarea de 


calcul: 
1 


-— 5 — 1,0001 - 107%. 
l— n 


Să presupunem acum că x, este aproximația de ordin n a rădăcinii simple, 
a ecuaţiei f(x) = 0. 
Atunci, din teorema lui Lagrange: 
Flota) = a — 0) 


unde !, aparţine intervalului generat de 4, și a. 


Dacă |f'()| > M, pentru x e J atunci: 


f(x 

| ie „ea DAE (14) 
M, 

Deci avem o estimare a erorii numai funcție de valorile lui /(x,) (ceea ce 

este foarte convenabil, de multe ori, dacă cunoaștem comportarea primei 


derivate). 


Să notăm cu Ț(x,) valoarea calculată a funcţiei f în x, și fie: 


(ac) — Plata) + dom), 30) să. 
Vom putea deci exprima cît de bine este aproximată soluția în funcție 


de 5. 
Situaţia cea mai bună ar fi să determinăm un x, astfel ca Îl) = 0. În 


acest caz valoarea exactă a funcției f(x) satisface inegalitatea: 


f(x] < 3 
Presupunînd că derivata întîi f'(4) nu variază foarte mult în apropierea 


soluției x = «, din (14) deducem: 


bă 3 
| — as — 8 ex unde e, = i (15) 
M, LF'(e)| 
Urmează că cea mai bună margine a erorii pentru orice metodă este e, 
care poartă numele de precizia ce se poate atinge pentru rădăcina a. 


53 


Dacă |/'(a)| este mic atunci e, este mare şi în acest caz convergența 
șirului 4x,) este slabă (spunem că problema este şău condiționată). 
Dacă privim figura III. 4, : 


E 


Fig. III. 4 


observăm că o mică schimbare în funcţie produce o mare variaţie în evaluarea 
rădăcinii. 

Pentru oprirea practică a procesului iterativ nu se folosește s, (nefiind 
cunoscut) ci următorul. criteriu: 


| Xa — Xa SS |ăa — Xacul (16) 
| E. Apar] Se d. 


$ 4. Metoda Newton-Raphson 


Metoda Newton-Raphson este o metodă de rezolvare a ecuaţiilor neli- 
niare bazată pe aproximarea curbei y = f(x) cu tangenta sa în punctul 
(20, f(%0)), apoi următorul punct 4, se obține intersectînd tangenta la grafic 
în X9 cu axa Ox. Procedeul continuă (fig. III. 5). Scriind ecuația tangentei în 
(o, f(0)) la graficul funcţiei y = f(x) obţinem: 

y — flo) = f'(0) (ăi — 0), de unde pentru y=—0 avegi (presupunînd 


Po) z 0): 
ao d Aa 

Deci: 
iu ma 0, | d 


Se formează. astfel şirul iteratelor xp, 2, Xp ses Xp. dat de algoritmul: 


ft), A 
f(x) 


Xa Ba hp, unde h, = — 


! i Fig. UL.5 


Funcţia /(x) se aproximează prin tangenta sa în punctul (x, /(4,)) şi 
Xa este abscisa punctului de intersecţie a tangentei cu axa Ox. 

Șirul (4) converge (nu dăm demonstraţia aici) către o rădăcină a ecuației 
/(%) = 0, dacă această rădăcină este presupusă simplă și alegerea inițială xo 
este suficient de apropiată de ea. 

Dăm mai jos un criteriu de convergență, ușor de aplicat. 

Teoremă. — Presupunem că /'(3) 40 şi /'"'(x) nu-și schimbă semnul 


în intervalul [&, 9] iar f(a) :/(5) <0. Dacă De] <b—a, Jb) 
) /'(b) 


'(a 
< b— d ae metoda Newton-Raphson cveldă (CA E orice alegere inițială 
x € la, 0]. « 

Pentru înțelegerea acestei teoreme a se vedea figura III.6. 

Exerciţiu. Verificaţi că teorema precedentă se poate aplica pentru rezolvarea 
ecuației x? — 1,5=—0 pe intervalul [1, 2]. 

Din şirul (x, dat de relaţia (19) noi putem calcula, evident numai un 
număr finit de termeni — de aceea trebuie să spunem cînd oprim procesul 
de calcul (19) astfel încît să avem o aproximaţie dorită a rădăcinii căutate. 


Criteriul de oprire (vezi: programul de calcul) va fi următorul: 


y | ai 
5 
AD 


Cînd 4, devine mai mic decît eroarea admisă în soluţie atunci procesul 


„de calcul se opreşte şi x, devine soluție aproximativă a ecuației f(x) = 0. . 


Dăm în figura III. 7 organigrama și programul FORTRAN pentru metoda 
Newton-Raphson: 


IN 


Tipăreşte : 


f(x.) 
fig i d Fat 


x 


vel 


Tipăreşte : 


Ni X ie a fiii) fi) 


Fig. III. 7. Organigrama pentru metoda Newton-Raphson 


Programul FORTRAN 


pentru metoda Newton- Raphson, 


READ (105, 1) XI, EPS 
FORMAT (GE 12.5) 
HF = F(XI) 
HI = HE/ED(XI) 
SCRI) = ST EEII 
IF(ABS(HI) — EPS) 10, 10, | 
il WRITE (108, 2) XI, XIL, HE, HI 
2 FORMAT (3X, GE 12.5) 
XI = XII 
GO TO 3 
10 HF = F(XII) 
WRITE (108, 2) XII, HE, HI 
STOP 
END 


am 
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FUNCTION F(X) 

F = ...aici se scrie programul pentru f(x) 
RETURN 

END 

FUNCTION FD(X) 

FD = „aici se scrie programul, pentru f(x) 
RETURN 

END 


Programul de mai sus a fost utilizat pentru rezolvarea ecuației: 
f(x) = sin x — (/2)2 = 0, 


S-a pornit cu aproximația inițială 4 = 1,5. Procesul iterativ a condus 
la soluție (cu precizia de 1 - 104) în numai patru pași. 

Exerciţiu, Rezolvaţi aceeași problemă pe calculator completind Programul 
Pentru funcţiile F(X) şi FD(X). Folosiţi în calcul 5 zecimale şi porniți cu apro- 
ximația xp = 1,5. 

Tapăviti la fiecare iterație următoarele valori: 


Nr. iteraţiei pă Xa Fl) hu | 


Exerciţii 

1. Pentru a intui convergența sau divergența unei metode iterative de 
rezolvare a ecuației + — F(x) faceţi reprezentări grafice pentru cazurile 

e 19 (0) e (Do 0/18 (0) ci 15 D003) Sale A) se la 

2, Explicaţi cu cuvintele dumneavoastră de ce putem spune că metoda 
Newton-Raphson este o metodă de „liniarizare“ plus iterare. - 

3. În cîţi pași se rezolvă ecuația x — 3 =0 prin metoda Newton-Raphson? 
Este importantă alegerea inițială ? 

4. Arătaţi că ecuaţia: 4 


ad — Sa? — 45 +50 = 0 are o rădăcină cuprinsă între 1 și 2. 


Calculaţi acea rădăcină prin metoda Newton-Raphson alegînd x = 1 (calcu- . 
lul cu 6 zecimale corecte). 


5, Calculaţi V2 prin metoda lui Newton-Raphson cu 6 zecimale exacte. 
6. Calculaţi 5 utilizînd metoda Newton-Raphson: 


i i 5 
Aaa 2%, AF 2 
3 Ev 
(calculele cu 6 zecimale corecte). 
7. Găsiţi cea mai mică rădăcină pozitivă a ecuaţiei: 
7 a ME dies 


(cu şase zecimale). 


8. În calculul rădăcinii «a ecuaţiei « = F(x) cu o eroare mai mică decît 
0,5. 104 s-a obținut: | 


pi 


Xa = 0,43789 și us — 0,43814. 


57 


„Se ştie că |F'(x)| s 0,4. Ce număr de iterații sînt necesare pentru a fi 
siguri că am atins acurateţea cerută? 
9. Scrieți un program de calcul (subrutină) pentru rezolvarea ecuației 
x — F(x), iterativ, %uu = Fa), cu xo aproximarea inițială și oprirea cal- 
culului dată de: . 


| Xa ua | 
10. Aplicați programul precedent la găsirea rădăcinilor ecuației sii și 
3x — cosy =0. = 


11.a) Rădăcina reală a ecuaţiei x? = x -| 4 poate fi scrisă: 


3 
2+ 321 |e- 521 


utilizați această expresie și calculați cu 4 zecimale corecte pe a. 
b) Calculaţi pe « cu metoda -Newton-Raphson tot cu 4 zecimale corecte, 
plecînd de la x — 2. | iu 
12. Calculați cu 4 zecimale. soluţia nenulă a ecuaţiei 


pr E. je (eiu- 
Exerciţii facultative (cap. LII) | 
__1) Fie g(x) o funcţie reală de variabilă reală și, xoe R. Vom spune că 
&(2) se comportă bine în xp şi vom scrie'g e FB(x0), dacă există, două inter- 
vale Î, = [%o, xo + e) și În = (o — e, X9], e > 0 astfel încît g() este convexă 
sau concavă pe Î, și independent convexă sau concavă pe Io. i 104 

a) Arătaţi că dacă g'(49) 2 0 atunci ge FB(axy). | 

b) Aparține funcţia e(x) = x? sin 1/x lui FB(0)? | 

c) Fie f: R-— R cu derivată continuă pe 1 = (30 — e, xo-l- e) pentru 
un anumit € arbitrar. Presupunem că: 


2(£) = Alto + 2) — f(ao — 2) are proprietatea ge FB(0). 
15.98 


d() lao + h) Pe —h) A 


arătați că pentru k > 0 suficient de mic 4(4) converge monoton către /'(%0). 


d) Pentru g(?) definit la punctul c) cu derivată de ordinul al doilea con- 
tinuă pe / = [0, e), e > 0 arbitrar și g, e” e FB(0) arătaţi că: 


F(h) = leo == sa 


este monoton crescătoare pentru / suficient de mic. 


Că 
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. 


e) În ipotezele de la punctul d) presupunînd în plus că avem un șir de 
numere pozitive monoton descrescător (4, care satistace pentru orice B > ii 
relația: i 

hi 


fag tt 0 Doo 
arătați că: 
pentru toți 2 > N (N un anume număr întreg) avem: 
Ida) — dh) | < dh) — dl. 
J) Tabelaţi derivatele funcțiilor elementare e, sin x, tg x, ja pe [0, 1] 


într-o sută de puncte echidistante utilizînd algoritmul dat de schema logică 
următoare (d(/) este dat la punctul €) schema III. 8. 


dfp)=1 200) Mhp 
M(h;)=-109 (544) 


pdreea 
, 
fa D(h,) 


Fig. 11.8 


x — 0 porniţi cu n = 0,01 +4. P. 
De asemenea am notat cu: 


DH) = (PA GEB) —FAXON)oO LO B) 


unde literele mari reprezintă reprezentarea în calcul a literelor mici iar ope- 
rațiile încercuite reprezintă aproximarea mașinii la operaţiile aritmetice, 


Precizăm că $ poate fi luat 4, P este precizia cu care lucraţi, iar dacă 


! 
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Răspunsuri 
” Capitolul I 
1. 8, 32768, 2660. 


2. a = 0,3333 - 10% eroarea =— — 1/3 + 1074, corecția = 1/3 - 1074. 
3. da. 
4. numai cînd în reprezentarea a = m - 10% -- m, - 107? (vezi $ 5) 
0 < |n] <0;5. 
5. 0,5: 10-6 pentru folosirea rotunjirii și ] - 10-% pentru trunchiere. 
6. 43, 56, 28, 39, 80. | 
Ze fial 1 doal ji 12097 pas be SA 25 42 ui 2A za al 2 200 2 e 8 IE 


1255; - 255 = 128, deci se adaugă 128 la exponentul real 


> 


q e [— 128, 127]. 
8. a) 0,2135 - 1051 - 0,3142 - 105: — 0,0671 - 1052 = 0,6710 - 10%), 
d) 0,2135 - 105% — 0,2120 - 1050 — 0,0015 - 10% = 0,1500 - 10%, 
e) 0,1214 - 1051 + 0,2224 - 1052 — 0,0121 - 1052 -+ 0,2221 - 10% = 
= 0,2342 - 1082, 
9: 0,087; 0,009; 0,0009 a treia. 
10. 0,1. 


11. ut (2 + 3) = 0,12115, 
(da -F Xa) -k a = 0,12114, 
(aa A 38) + Xa = 0,1214. 


n n— 
12. a) Avem vw, Spa = ( ai PE, da == /ns 
- | moda) 
b) vo = In 6 — In 5 = 0,182, 
1 = 0,0090, 
Va = 0,050, 
Va = 0,083, e 
Va = — 0,165 (absurd). 


c) Eroarea de rotunjire 5 - 10-4 se amplifică la fiecare pas astfel că 
obţinem în y, o eroare de 625 - 5 : 104 = 0,3125. 


13. yp = 0,019, ys = 0,028, ya = 0,058, 
Yp = 0,021; ya = 0,034, m = 0,088, 
6 = 0,025, yg = 0,043, yo = 0,182 (corect). 


Deci avem stabilitate numerică. 


Capitolul II 
1. Se exprimă ctg « în funcție de parametrii directori ai celor două 
drepte. 


2,. Se foloseşte expresia numărului de condiţionare dat de lungimea 
euclidiană a unei matrice. 


3. cond (4) = 10.967. Soluţia sistemului inițial « = 1; y = 1. Soluţia 
sistemului perturbat % — —12,9; 3 = —20. Perturbarea soluției x — % = 
— 139; y—y=a2l. : 


4. Nu este nevoie. 
5. x — 6,154047, x — 4,313317, Xa = 3,240207. 


6. După cinci iterații obținem: 


x = 0,49854, 
Xp = 0,74927, 
sp = 0,24921, 
sa — 0,49963. 
7. conda (4) > 15:c. 
8. Soluția sistemului inițial z= — 1; y=0. Soluţia sistemului per- 


turbat 7 = 1; = 
- | —2 
Se notează A Sa 
, —2 4,01 : 
se înmulțește cu 0,01 și se deduce că apriori putem avea o modificare de 
0,01 în maximum pe componente. 


şi se calculează ||A-!|„ care număr 


Exerciţii facultative. b) Se obține relația de recurenţă 


Â, = 2 cos 6 A Ana 


Capitolul III 

3. Un singur pas. Nu. 

4. fa) — m — 5x2 — 45x + 50. Se verifică relaţia: 
JUD - 12) <0. 


5. Aplicind metoda Newton-Raphson obținem: x, = - (2/0 Pe a). 


7. % = 4,49341. 
8. Se aplică rezultatele din $ 3. 


9. SUBROUTINE ITER (FUNC, X0, EPS) 
X1 = PUNC(X0) 
X2 = FUNC(X |) 


6. X3 = FUNC(X2) 

SEM 2 (3 2 [020 Xa) 

A = XM.(X3—X2)/(1— XM) 

IF(A.LTEPS) GO TO7  * 

cp ai) 

2 ==,X3 

GO TO 6 

WRITE (108,8) X3 <a pata 

FORMAT (IX, E 12.5) 

4 STOP j die 
END? : LE 


a 


a ș cos 4, 4 : PIE rr 
10. Luaţi 40 == 0,317 şi su = * veţi obține rădăcina cuprinsă între 


0,316749 şi 0,316754. 


11. 1,7963. 


12. Considerăm. f(x) — x — 1 he =—0, 
Metoda Newton ne dă: aa 


Aa Xa A LE ha Ş (37 l i ez) /(1 "Fi Ze); 


Pentru -xp = 0,8 obţinem x, = 0,79682 deci hy = — 105 deci rădăcina 
este e — 0,7968. | 


Răspunsuri la exercitiile facullative (cap. II) Si 
c) avem: d(h) = e, prin derivare obținem: 
va ! 1 2h 4 : 
; "(h) - h— gh 2 (h) - h — (e(h) — e(0) 
TIPA pe A : EI) eh (e) 240) 
2 - hi „p2h 


_ gi) —e 
a hi 
Deoarece g e FB(0), pentru / suficient de mic (d(h)) are exact un singur 
semn apoi din fi > fz(/n și fn» fac parte din acei / pentru care raționamentul 
de mai sus merge) deducem: 


(ha) — d(n) = ii d(s)ăs 


„0 <1<h (am aplicat formula lui. Lagrange pe LO, h]). 


i 


.. 


de unde rezultă monotonia, apoi 


lim d(4) = f'(%0) pentru / cu derivată continuă, 
h->0 ) 


d) Deoarece g e FB(0) putem presupune că g este convexă pentru k 
suficient de mic, atunci g' este crescătoare pentru / suficient de mic, deci: 


Fi) = e. 


h 
Prin derivare obținem: z 
F*(h) a e"'(h) e (2) e 2(P)/h Ş 2"(h) 8 es 5. O<i<h 
, Aa 
= g"(h) — e" (SU — 9h), t< sh. zi 


Deoarece g' e FB(0) urmează că g'” este monotonă. Cum g este presupus 


convex urmează că g” > 0 iar g(/) = f(ao 4-1) — /(Xo — 1) ne dă g”(0) =0. 


Dacă g” este monoton crescătoare atunci F'(h) > 0 în timp ce e” este 


„monoton descrescătoare urmează că F'(/) = 0 căci g'(0) = 0 şi e” >0. 


e) Urmînd punctul c) putem scrie: 


la) = da =63 


ii 2 h 


FT pppeca st. | dh 


Ne], 


conform lui d) pentru 7 3 N, N un anume întreg pozitiv, avem: 


at) — aţa) > IE El = pi aa = 


hi 


] Ed In (fesa/h) fi lg'(h,) = le 


Similar avem: 
LR , 
la) = Aha) -4 le. pala as 
LI ze [A 


IE) — gt)! 


< în (hihi 


de unde rezultă concluzia cerută. 
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